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ABSTRACT 

Timing observations from the Parkes 64-m radio telescope for 165 pulsars between 
1990 and 2011 have been searched for period glitches. Data spans for each pulsar 
ranged between 5.3 years and 20.8 years. From the total of 1911 years of pulsar rota- 
tional history, 107 glitches were identified in 36 pulsars. Out of these glitches, 61 have 
previously been reported whereas 46 are new discoveries. Glitch parameters, both for 
the previously known and the new glitch detections, were measured by fitting the tim- 
ing residual data. Observed relative glitch sizes A^ g /V range between 10 -10 and 10 -5 , 
where v = l/P is the pulse frequency. We confirm that the distribution of Av g /i> is 
bimodal with peaks at approximately 10~ 9 and 10~ 6 . Glitches are mostly observed in 
pulsars with characteristic ages between 10 3 and 10 5 years, with large glitches mostly 
occurring in the younger pulsars. Exponential post-glitch recoveries were observed for 
27 large glitches in 18 pulsars. The fraction Q of the glitch that recovers exponentially 
also has a bimodal distribution. Large glitches generally have low Q, typically just a 
few per cent, but large Q values are observed in both large and small glitches. Ob- 
served time constants for exponential recoveries ranged between 10 and 300 days with 
some tendency for longer timescales in older pulsars. Shorter timcscalc recoveries may 
exist but were not revealed by our data which typically have observation intervals of 
2-4 weeks. For most of the 36 pulsars with observed glitches, there is a persistent 
linear increase in i> (i.e., decrease in the slow-down rate \i>\) in the inter-glitch interval. 
Where an exponential recovery is also observed, the effects of this are superimposed 
on the linear increase in v. In some but not all cases, the slope of the linear recovery 
changes at the time of a glitch. The v values characterising the linear changes in v 
are almost always positive and, after subtracting the magnetospheric component of 
the braking, are approximately proportional to the ratio of |z>| and the inter-glitch 
interval, as predicted by vortex-creep models. 
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1 INTRODUCTION 

Pulsars are thought to be highly-magnetised, rapidly- 
rotating neutron stars. They are remarkably stable ro- 
of general relativity 



tators, which has enabled tests 
■ ' it —\ 

(Kra mer et al.l|200(: ). searches for gravitational waves (e.g 
lYard 



Yardlev et al.l 12011 ) and the establishment of a pulsar 



timescale (|Hobbs et al.ll201ll ). These results have been be- 
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ing obtained by the technique known as "pulsar timing". 
The pulsar timing technique allows observed pulse times-of- 
arrival (ToAs) to be compared with predicted arrival times. 
The predicted arrival times are determined using a model of 
the pulsar's rotation, position, orbit etc. The differences be- 
tween the actual and predicted pulse arrival times are known 
as "timing residuals". Timing residuals can be induced by 
an inaccuracy or omission in the parameters in the timing 
model or by the timing model not including all phenom- 
ena affecting the propagation of a pulse from the pulsar to 
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the observer. The timing residuals for some pulsars are very 
small. For instance, PSR J0437— 4715 has an rms residual of 
75 ns over several years llManchester et al.ll2 012). However, 
most pulsars are not so stable. iHobbs et al.l (2010) analysed 
the timing residuals of 366 normal and recycled pulsars on 
timescales longer than 10 yr. They found that in most cases 
the residuals comprise low- frequency structures. For young 
pulsars, the timing residuals were further found to be dom- 
inated by recovery processes from glitch events. 

A glitch is an abrupt increase in the pulse frequency 
v = 1/P of a pulsar, often followed by an exponential re- 
covery t oward the extrapo lation of the pre-glitch pulse fre- 
quency (|Bavm et al.ll 1963 ). Post-glitch behaviours generally 
exhibit another recovery process which is characterised by 
a linear increase in i> or decrease in slow-down rate \i>\. 
This often extends from the end of the initial exponential 
recovery until the next glitch event. Such "linear -decay" 
processes were first o bserved in the Vela pulsar (IPownsI 
ll98ll ; lLvne et al.lll996h, and were subsequently seen in other 



( Yuan et al. mOlOal ). The first known glitch 



detected in the Vela pulsar jRadhakrishnan fc Manchester! 
Il969l : iReichlev fc Downs! ll969l V Since then more than 350 
glitch events have been observed in about 120 pulsars. Glitch 
databases a re now available: the A TNF Pulsar Catalogue 
glitch table dManchester et al . 2005 f] and the Jodrell Bank 
Glitch Catalogue ( Espinoza et al. " 1201 lah PI Since the orig- 
inal glitch discovery, the Vela pulsar has been observed 
to undergo 15 further glitch events, most of which have 
Au g /v ~ 10 -6 . In contrast, the Crab pulsar has been ob- 
served to have Av g /u ~ 10~ r — 10~ 9 for most of its glitch 
events. Most glitches have been observed in relatively young 
radio pulsars but they have also bee n observed in magnetars 
l|Woods et al.ll2004l; iDib et al1l2007l) a nd even in a millisec- 
ond pulsar (jCognard fc Backed 120041 ). Observed fractional 
glitch sizes range from ~ 10~ 10 to ~ 10 -5 , but it is impor- 
tant to note that the low end of this distribution is strongly 
limited by observational selection. 

The increase of the pulse frequency during a glitch 
is usually unresolvable and exponential recover ies typically 
have timescales of ten to a few hundred days (jWang et al.l 
l2000l ; lYuan et al. I l2010al ). However, intensive observations 
of glitch events in the Crab and Vela pulsars have shown 
that 1) the rising edge of the pulse frequency can some- 
times be resolved into multiple components (|Lvne et al.l 
1 1992ft and 2) very short exponential decays can occur 



( Dodson et al.ll2002ft . For two pulsars, sinusoidal oscillations 
have been observed in timing residuals a fter glitch events 
(|McCulloch et al.lll990l ; lYuan et alfolObl ). 

"Slow glitch" events ha ve been observ e d in PSR 
B1822-09 and other pulsars ilZou et all |2004| ; IShabanoval 
l2005l . l2007l : lYuan et all l2010al ). Unlike normal glitches, a 
slow glitch builds up over several hundred days, and the 
increased pulse frequency is usually maintained until the 
next event. This corresponds to a fluctuation in charac- 
terised by_jui_jrrn2ulsiyedecre^ in- 
crease. IHobbs et all (|201Cft and lLvne et~al1 (|2010h suggested 
that slow glitches are a manifestation of the il v switching" 
observed in some pulsars. 



1 http: / / www.atnf.csiro.au/research / pulsar / psrcat / glitchTbl.html 

2 http:/ /www. jb.man.ac.uk/pulsar/glitches. html 
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Figure 1. Period-period-derivative (P - P) diagram showing the 
pulsars in our sample where no glitch is detected (x), a glitch is 
detected (A) and a glitch was detected prior to our observations 
(O)- Other pulsars are marked with a dot. Data are from the 
ATNF Pulsar Catalogue. 



Glitches are thought to be 



either by 



tri ggered 

the neutron-star cru stquakes (e.g., iRudermanl Il99ll : 
iRuderman et al] 1 19981 ) or by the sudden transfer of an- 
gular momentum from the faster-rotatin g crustal neutron 
superfluid to the rest of the star (e.g., lAnderson fc Itohl 
1 19751 ; IRudermanl Il976l : lAlpar et al.lll98ll ). The post-event 
exponential recoveries have been explained as the re- 
establishment of an equilibrium between pinning and un- 
pinni ng in a vortex-creep region inter ior to a neutron 
star (jAlpar et al.lll993l : iLvne et al.ll2000l ). Fractional glitch 
sizes Avgjv show a bimodal distribut i on with peaks at 
~ 10~ 9 and ~ I P" 6 l|Lvne et all 1200(1 IWang et all l2000l : 
Yuan et al . 2010a). Using a sample containing 315 glitches, 
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2011aft confirmed this bimodal distribution 
and also found that the rate of glitch occurrence peaks for 
pulsars with a characteristic age (r c = P/(2P) of about 
lOkyr. They also showed that, on average, nearly one per 
cent of the spin-down is reversed by glitches for those pulsars 
with a slow-down rate \i>\ between 10~ 14 to 10" 11 s~ 2 . 

Even though glitch events and their subsequent recov- 
eries have been extensively studied, theoretical predictions 
have been unable to model fully the timing residuals induced 
by a glitch event. Theoretical models also cannot yet explain 
why some pulsars exhibit a large number of glitch events, 
whereas other pulsars wi th similar characterist ics have never 
been observed to glitch. iMelatos et al.l (|2008l ) showed that 
the waiting-time sequences of the glitches in seven pulsars 
followed a constant-rate Poisson process, which suggests 
that a neutron star could be a system fluctuating around 
a self-organised critical state. 

For this paper, we searched a total of 1911 yr of pulsar 
rotational history for glitch events. In §2, we describe our 
observations. In §3, we present our method for determining 
glitch parameters. Our results are shown in §4, and discussed 
in §5. We conclude the paper in §6. 
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Figure 2. Illustration of the identification of glitch events with 
two typical sizes 10 — 6 and 10 — 9 . Sub-plot a) shows the effect on 
the timing residuals of the large glitch that occurred in the Vela 
pulsar around MJD 53193. The pre-glitch solution contains v, v 
and v and phase-coherency is broken by the glitch. Sub-plot b) 
shows timing residuals relative to a pre-glitch solution containing 
v and i> for a small glitch in PSR J0834— 4159 that occurred at 
MJD ~ 53415. In this case, phase coherency is maintained after 
the glitch although there are several phase wraps. In each plot 
the vertical dashed line indicates the glitch epoch (from Dodson 
et al. 2004 for the Vela glitch). 



2 OBSERVATIONS 

The observations of 165 pulsars analysed in this paper 
were obtained using the Parkes 64-m radio telescope be- 
tween 1990 January and 2011 January. Almost all obser- 
vations were at radio frequencies near 1400 MHz, in the 
20-cm band. For 1990 - 1994 the "H-OH" receiver was 
used with an analogue filt erbank having 64 x 5 MHz chan- 
nels for each polarisation (| Johnston et al.l ll992). The data 
were summed and high-pass filtered before one-bit digitisa- 
tion. From 1994 - 2001 most data were obtained using the 
"FPTM" digital filterbank (jSandhu et all 1 19971 ) with one 
(later two) 128 MHz bands for each polarisation. Up until 
1997 the H-OH receiver was used. From 1997 on, most ob- 
servatio ns used the centre beam o f the 20-cm Multibeam 
receiver (|Stavelev-Smith et al.l[l99o1 ) although a few were at 
higher frequencies between 1700 and 2200 MHz using other 
receivers. Between 1997 and 2007 observations were usu- 
ally made with the analogue filterb ank system used for the 
Parkes Multibeam Pulsar Survey (|Manchester et al]|200ll ) 
which has 96 x 3 MHz channels for each polarisation. Inter- 
vals between observing sessions were typically 2-4 weeks 
and a ToA was obtained for most of the pulsars each ses- 
sion. Some pulsars with lower priority were observed less fre- 



quently. Observation times per ToA were normally between 
one and ten minutes. 

From 1991 to 2000 these observations were in part used 
to support the Energetic Gamma-Ray Experiment Tele- 
scope ( EGRET) aboard the Compton Gamma-Ray Obser- 
vatory (|Thompsonll200St ). Between 2007 and 2011, obser- 
vations were obtained with the primary goal o f supporting 
the Fermi Gamma-ray Space Telescope mission IjSmith et al.l 
120081 ; IWeltevrede et al.ll2010l ). Observing sessions are sepa- 
rated by approximately four weeks, each lasts for 24 h, al- 
lowing ~170 pulsars to be observed. The centre beam of the 
Multibeam receiver is used at 1369 MHz with a bandwidth 
of 25 6 MHz. Digital filterbank systems (|Manchester et all 
l2012h were used to record the data, with integration times of 
2 to 20 min for each pulsar to ensure a signal-to- noise ratio 
larger than five. A few ob servations were tak en as part of the 
PULSE@Parkes project (|Hobbs et al.ll2009l ). These data are 
available for downl oad from the Parkes pulsar data archival 
jHobbs et al.ll201ll ). 

In Table [TJ we summarise the properties of the 165 pul- 
sars. The pulsar name, pulse period P, period derivative P, 
dispersion measure (DM), data span in Modified Julian Day 
(MJD) and years, and number of observations are presented. 
The final column in this table indicates whether the pulsar 
has never been observed to glitch (N), has been detected 
to glitch during our observations (Y) or has been reported 
to glitch prior to our observations (P). Figure [T] shows the 
pulsars in our sample on the period-period-derivative (P - 
P) diagram The P and P data are from the ATNF Pulsar 
Catalogu <Q ^Manchester et alj|2005l ). The identification and 
modelling of glitch events will be introduced in detail in the 
following section. 



3 DATA ANALYSIS 

Off-line data re duction used the psrchive pulsar data 
analysis system (| Hot an e t al. 2004). Each observation was 
summed in time, frequency and polarisation to form a to- 
tal intensity pulse profile. In order to determine the pulse 
time-of-arrival (ToA), each of the total intensity profiles 
was cross-correlated with a high signal-to-noise ratio "stan- 
dard" profile. Timing residuals were for med using the pul- 
sar timing software package Tempo2 (|Hobbs et al.l 120061 ; 
lEdwards et al.ll2006l ). with the Jet Propulsion Labo ratories 
(JPL) planetary ephemeris DE405 |Standish|[l998l ) to cor- 
rect the local ToAs to the solar-system barycentre. Each ob- 
served ToA was first referred to terrestrial time as realised by 
International Atomic Time and subsequently to Barycentric 
Coordinate Time. For each pulsar, Tempo2 was used to find 
a set of parameters that provided a phase-connected timing 
solution. The solution contains the pulse frequency v and its 
first derivative i>. The pulse frequency second derivative was 
only fitted when a cubic structure in timing residuals could 
be seen after fitting for v and v. In some cases, particularly 
for pulsars that had glitched or have large amounts of tim- 
ing noise, it was not possible to obtain a phase-connected 



3 http://data.csiro.au 

4 http://www.atnf.csiro.au/research/pulsar/psrcat/; catalogue 
version 1.43 
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Table 1. Pulsars in our sample observed at the Parkes Observatory with a data span larger than 5 yr. 



PSR ,1 PSR BP P DM Data range Data span No. of ToAs Glitched? 
(b) (1Q- 15 ) (cm- 3 pc) (M.ID) (yr) (Y/N/P) 1 



J0108 — 1431 









807565 





08 


2 


.38 


49373 




55144 


15 


-8 


135 


N 


,J040 1 — 7608 


B0403 


— 76 





.545253 


1 


54 


21 


,60 


53033 




55144 


5 


.8 


58 


N 


,10536 — 7543 


B0,)3H 


— 75 




.245856 





56 


17 


,50 


48957 






16 


.9 


103 




J0630 — 2834 


B0628 


— 28 




. 244419 


7 




34 


,47 






55144 


9 


.9 


62 




J0729-1448 









.251659 


113 


29 


92 


,30 


51896 


_ 


55429 


9 


7 


89 


Y 


J0738 — 4042 


B0736 


-40 





.374920 


1 


62 






52995 




55144 


5 


9 


73 




J0742 — 2822 


B0740 


-28 





. 1 C>(>7(>2 


16 


82 


73 


.78 


49364 




..)..)..> <9 


17 





481 




J0R34 — 4159 









.121116 


4 


44 


. i ■ 


.50 


i . 




55145 


10 




120 




J0835 — 3707 









. 5 '114(14 


9 


.78 




,30 


509 10 




lj3948 


8 


2 


62 




J0835-4510 


B0833 


-45 





.089328 


125 


.01 


67 


,99 


49608 


Z 


55172 


15 


2 


667 


Y 


,10855 — 4644 









. 064686 


7 


26 


238 


.20 


51158 




55144 


10 


.9 


225 


N 


.108") 7 — 4424 









.326774 




.34 


181 


,43 


51899 




55144 


8 


.9 


94 




,70901 — 4624 









.441 9 9 5 


87 


49 


198 


,80 


5 (184 9 




55144 


11 


.7 


127 


N 


,10905 — 5127 









.346287 


24 


90 


196 


,43 


49363 




55145 


15 


.8 


101 


Y 


J0908-4913 


B0906 


-49 





.106755 


15 


15 


180 


.37 


48957 




55182 


17 


.0 


271 


N 


J0940 — 5428 









.087545 


32 


-87 


134 


.50 


509 41 




55144 


11 


.5 


163 




J0942 — 5552 


Bil'.l 111 


-55 





. (>(> 4307 


22 


85 


180 


.20 


48928 




53948 


13 


7 


181 


N 


,10954 — 5430 









.472834 


43 


.91 


200 


,30 


5094(1 




55182 


11 


.6 


106 




J 1012 — 5857 


B1011 


— 58 





.819911 


17 


69 


383 


.90 


5 5 3 6 




53948 


9 


.3 


64 




J1015-5719 









.139882 


57 


37 


278 


.70 


51215 


Z 


55182 


10 


.9 


166 


N 


,J 1016 — 5819 









.087834 





.70 


-' 


,10 


509 4(1 




55144 


11 


.5 


77 




,J 1016 — 5857 









. 107386 


80 


83 


394 


.20 


5 1 2 9 9 




55429 


11 


.3 


250 




J 1019 — 5749 









. 1 62499 


20 


08 


1039 


.40 


51158 




55182 


11 


.0 


85 




J 1020 — 6026 









. 140480 


6 


74 




,00 


52854 




55182 


6 




67 




J1038-5831 


B1036 


-58 





.661992 


1 


25 


72 


.74 


50536 




53948 


9 


.3 


61 


N 


J 1043 — 6116 









.288602 


10 


40 


449 


,20 


51158 




55182 


1 1 


.0 


75 


N 


,1 1047— 6709 









.198451 


1 


69 


116 


,16 


50538 




53948 


9 


3 


62 


N 


,11048 — 5832 


B 1046 


— 58 





.123671 


96 


32 


129 


,10 


47910 




55183 


19 


.9 


353 


Y 


,1 1052 — 5954 









.180592 


19 


98 




.00 






55460 






92 




J1057-5226 


B1055 


-52 





.197108 


5 


83 


30 


,10 


49363 


Z 


55182 


15 


.9 


291 


N 


J 1 105 — 6107 









.063193 


15 


.83 


271 


,01 


49589 




55461 


16 


1 


297 


Y 


J 1 1 12 — 6103 









.064962 


31 


46 


599 


,10 


5 (185(1 




55207 


11 


.9 


178 




,71114 — 610(1 




-60 




.88082(1 










50538 




53948 










J1115 — 6052 












23 


228 








d20o 


11 


9 


97 




,11119-6127 









.407963 


4020 


22 


707 


,40 


50852 


Z 


55576 


12 


.9 


348 


Y 


J 1 123 — 6259 









.271434 


5 


25 


223 


,26 


5040(1 




55205 


13 


1 


182 


P 


,1 1 136 — 5525 


B1133 


— 55 





.364706 


8 


-- 


85 


,50 


5184 4 




53948 


5 


.8 


33 


N 


J 1 138 — 6207 









.117564 


12 


.48 


519 


,80 


5 (18 4 9 




55205 


11 


.9 


129 




,1 1 152 — 601 2 









.376570 


6 


68 


74 


,00 






53948 


7 


.5 


54 




,11156-5707 









.288409 


26 


45 


243 


.50 


51944 


z 


55205 


8 


.9 


67 


N 


,11216 — 6223 









.374047 


16 


.82 


786 


,60 


50851 




5 5 205 


11 


.9 


61 




J 1224 — 6407 


B1221 


— 63 





.216476 


4 


95 


97 


.47 


4 8 3 3 ( 1 




55205 


18 


.8 


661 




,1 1248 — 6344 









.198335 


16 


92 


433 


,30 


51260 




55205 


10 


.8 


78 




J 1301 — 6305 









.18 4 5 2 8 


266 


75 


374 


.00 


50941 




55104 






177 


Y 


J1305-6203 









.427762 


32 


14 


470 


,00 


50940 




55205 


11 


.7 


58 


N 


,11316 — 6232 









. 342825 


5 


30 


983 


.30 


49589 




53948 


11 


9 


163 




J 1 320 — 5359 


B1317 


— 53 





.279729 


9 


25 


97 


.60 


50530 




55205 


12 


.8 


178 




,1 1327— 6400 









.280678 


31 


18 


oso 


.90 


50940 




55205 


11 


.7 


63 




,1 1328 — 4357 


B1325 


— 43 





.532699 


3 


.01 




,00 


50738 




53948 


8 


.8 


88 


P 


,11341-6220 


B1338 


-62 





.193340 


253 


11 


717 


.30 


49540 


z 


55461 


16 


2 


265 


Y 


,11349 — 6130 









. 2 5 9 3 6 3 





12 


. '-4 


,60 


509 40 




55205 


11 


7 


71 


N 


,1 1359 — 6038 


B 1356 


— 60 





.127501 


6 


.34 


293 


,71 


48330 




55205 


18 


.8 


661 




,1 1412 — 6145 









.315225 


98 


66 


514 


. 70 


50850 




55461 


12 


6 


159 


Y 


,1 1413 — 6141 









.285625 


333 


.44 


677 


,00 


50850 




55461 


12 


.6 


198 


Y 


J1420-6048 









.068180 


S3 


17 


358 


.80 


51100 




55461 


11 


.9 


272 


Y 


,1 1452 — 585 1 









.386625 


50 


.71 


262 


.40 


51088 




55205 


11 


.3 


120 




.J 1452 — 6036 









.154991 


1 


45 


349 


.70 


51302 




55461 


11 


4 


93 


Y 


,1 1453 — 6413 


B1449 


— 64 





.179485 


2 


-75 


71 


,07 


5 (Hi 09 




55 205 


12 


4 


143 




,1 1456 — 6843 


B1451 


— 68 





.263377 





10 


8 


.60 


48330 




55205 


18 


.8 






J1509-5850 









.088922 


9 


.17 


140 


,60 


51214 


z 


55205 


10 


.9 


155 


N 


,11512 — 5759 


B 1508 


— 57 





. 128694 


6 


-85 




,70 


51527 




55205 


10 


1 


85 




J 1 5 1 3 — 5 9 8 


B 1509 


-58 





.150658 


1 5 3 6 


53 


2 5 2 


.50 


47913 




55205 


20 





3 8 4 


N 


J 1514 — 5925 









. 148796 


2 


-88 


194 


. 10 


51220 




55205 


10 


.9 


66 


N 


J 1515 — 5720 









.286646 


6 


10 


482 


,00 


51391 




55205 


10 


4 


54 


N 


J — ODiO 








78 1 9 
























J1524-5706 






1 


116049 


356 


.47 


833 


.00 


51101 




55205 


11 


.2 


110 


N 


J1530-5327 









.278957 


4 


68 


49 


,60 


51013 




55205 


11 


.5 


113 


N 


,11531-5610 









.084202 


13 


74 


110 


,90 


51215 




55461 


11 


6 


161 


Y 


J1538-5551 









.104675 


3 


21 


603 


,00 


51300 




55205 


10 


.7 


81 


N 


J1539-5626 


B1535 


-56 





.243392 


4 


85 


175 


.88 


49358 




55205 


16 





261 


P 


,11541-5535 









.295838 


75 


02 


428 


.00 


51300 




55205 


10 


.7 


69 


N 


J1543-5459 









.377119 


52 


02 


345 


,70 


50941 




55205 


11 


.7 


113 


N 


J1548-5607 









.170934 


10 


74 


315 


.50 


50941 




55205 


11 


.7 


127 


N 


J1549-4848 









.288347 


14 


11 


55 


.98 


49358 




55205 


16 


.0 


239 


N 


J1551-5310 









.453394 


195 


.13 


493 


.00 


51099 




55205 


11 


.2 


152 


N 


,11557-4258 









.329187 





33 


144 


,50 


50538 




53948 


9 


.3 


54 


N 


J1559-5545 


B1555 


-55 





.957242 


20 


48 


212 


.90 


49359 




53948 


12 


.6 


117 


N 


J1600-5044 


B1557 


-50 





.192601 


5 


06 


260 


,56 


50618 




55205 


12 


.6 


128 


N 


J1601-5335 









.288457 


62 


.37 


194 


,60 


50941 




55205 


11 


7 


118 


N 


J1602-5100 


B1558 


-50 





.864227 


69 


58 


170 


,93 


47913 




55205 


20 


.0 


246 


N 


,11611-5209 


B1607 


-52 





.182492 


5 


.17 


127 


.57 


51526 




55205 


10 


.1 


80 


N 


J1614-5048 


B1610 


-50 





.231694 


494 


94 


582 


.80 


47910 




55461 


20 


7 


413 


Y 


J1623-4949 









.725732 


42 


09 


183 


,30 


50851 




53975 


8 


.5 


57 


N 


J1626-4807 









.293928 


17 


48 


817 


,00 


50941 




55205 


11 


.7 


71 


N 


J1627-4706 









.140746 


1 


.73 


456 


.10 


52807 




55205 


6 


6 


84 


N 
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Table 1. — continued 



PSR J PSR BP P DM Data range Data span No. of ToAs Glitchcd? 
(s) (If)- 15 ) (cm- 3 pc) (M.ID) (yr) (Y/N/P) 



,11632-4757 









228564 


15 


07 


578 


00 


51216 




55205 


10 


.9 


92 


N 


J1632 — 4818 









813453 


650 


42 


758 


00 


50852 


_ 


55182 


11 


.9 


135 


N 


J1637-4553 


B1634- 


45 





118771 


3 


19 


193 


23 


50669 




55205 


12 


,4 


132 


N 


J1637 — 4642 









154027 


59 


20 


417 


00 


51393 




55205 


10 


,4 


86 


N 


J1638 — 4417 









117802 


1. 


61 


436 


00 


51633 




55205 


g 


.8 


82 


N 


J1638-4608 









278137 


51 


50 


424 


30 


51089 




55205 


11 


3 


75 


N 


J1640 — 4715 


B1636- 


47 





517405 


42 


03 


591 


70 


51528 




55205 


10 


,1 


39 


N 


J1643 — 4505 









237383 


31 


83 


484 


00 


52738 




55205 


6 


.8 


44 


N 


J1644 — 4559 


B1641- 


45 





455060 


20 


09 


478 


80 


47913 




55101 


19 


,7 


298 


p 


J1646 — 4346 


B1643- 


43 





231603 


112 


75 


490 


.40 


47913 




55273 


20 


.2 


305 


Y 


J1648-4611 









164950 


23 


75 


392 


90 


51216 




55205 


10 


.9 


73 


N 


J1649 — 4653 









557019 


49 


74 


332 


00 


51089 




55205 


11 


.3 


93 


N 


J1650 — 4502 









380870 


16 


06 


319 


70 


50941 




55205 


11 


.7 


82 


N 


J1650 — 4921 









156399 


1 


82 


229 


90 


52983 




55205 


6 


,1 


69 


N 


J1702 — 4128 









182136 


52 


34 


367 


10 


51089 




55205 


11 


.3 


87 


N 


.T1702-4310 









240524 


223 


78 


377 


00 


51223 




55461 


11 


.6 


125 


Y 


J1705— 1906 


B1702- 


19 





298987 


4. 


14 


22 


-91 


51901 




55206 


9 


.0 


75 


N 


J1705 — 3950 









318941 


60 


60 


207 


10 


51217 




55205 


10 


.9 


63 


N 


J1709 — 4429 


B1706- 


44 





102459 


92 


98 


75- 


69 


479IO 




55507 


20 


.8 


395 


Y 


J1713 — 3949 









392451 






342 


00 


51557 




54504 


8 


.1 


84 


N 


.T1715-3903 


* 







278481 


37 


69 


313 


10 


51217 




55205 


10 


.9 


111 


N 


J1718 — 3718 






3 


378574 


1613. 


59 


371 


10 


51244 




54859 


9 


.9 


100 


N 


J1718 — 3825 









074670 


13 


22 


247 


40 


50878 




55507 


12 


.7 


164 


Y 


J1721 — 3532 


B1718- 


35 





280424 


25 


19 


496 


00 


51879 




55205 


9 


.1 


96 


N 


J1722 — 3712 


B1719- 


37 





236173 


10 


85 


99 


50 


49363 




55205 


16 


.0 


197 


N 


.11723-3659 









202722 


8 


01 


254 


.20 


50851 




55205 


11 


.9 


102 


N 


J1726 — 3530 






1. 


110132 


1216 


75 


727 


00 


50681 




55205 


12 


.4 


147 


N 


J1730 — 3350 


B1727- 


33 





139460 


84 


83 


259 


00 


50539 




55507 


13 


.6 


182 


Y 


J1731 — 4744 


B1727- 


47 





829829 


163 


63 


123 


33 


48184 




55507 


20 


.0 


228 


Y 


J1733 — 3716 


B1730- 


37 





337586 


15 


05 


153 


50 


51893 




55205 


9 


.1 


80 


N 


,11734-3333 






1 


169008 


2278. 


98 


578 


.00 


50686 




55205 


12 


.4 


136 


N 


J1735 — 3258 









350963 


26 


08 


754 


00 


51393 




55205 


10 


.4 


64 


N 


J1737 — 3137 









450432 


138. 


76 


488 


20 


51157 




55507 


11 


.9 


83 


Y 


J1737 — 3555 


B1734- 


35 





397585 


6 


12 


89 


41 


52003 




53948 


5 


.3 


21 


N 


J1738 — 2955 









443398 


81 


86 


223 


40 


51158 




55205 


11 


.1 


63 


N 


J1739-2903 


B1736- 


29 





322882 


7. 


88 


138 


56 


50739 




55205 


12 


.2 


156 


P 


J1739 — 3023 









114368 


1 1 


40 


170 


00 


51879 




55205 


9 


.1 


98 


N 


J1740 — 3015 


B1737- 


30 





606887 


466 


12 


152 


.15 


50669 




55507 


13.2 


190 


Y 


J1745 — 3040 


B1742- 


30 





367429 


10 


67 


88 


37 


51901 




55205 


9 


.0 


121 


N 


J1752 — 2806 


B1749- 


28 





562558 


8 


13 


50 


37 


47911 




55083 


19 


.6 


177 


N 


J1756-2225 









404980 


52 


69 


326 


00 


51217 




54564 


9 


.2 


47 


N 


J1757 — 2421 


B1754- 


24 





234101 


12 


92 


179 


45 


51529 




55205 


10 


.1 


96 


N 


J1759 — 2205 


B1756- 


22 





460974 


10 


87 


177 


16 


51529 




53975 


6 


.7 


28 


N 


J1801 — 2154 









375297 


16 


00 


387 


90 


51218 




55205 


10 


.9 


61 


N 


J1801 — 2304 


B1758- 


23 





415827 


112 


93 


1073 


90 


47911 




55507 


20 


.8 


411 


Y 


.T1801-2451 


B1757- 


24 





124924 


127 


91 


289 


00 


48957 




55507 


17 


.9 


331 


Y 


J1803 — 2137 


B1800- 


21 





133667 


134 


36 


233 


99 


50669 




55530 


13 


3 


182 


Y 


J1806 — 2125 









481789 


121 


40 


750 


40 


51155 




55206 


11 


.1 


74 


p 


J1809- 1917 









082747 


25 


54 


197 


10 


50782 




55530 


13 


.0 


134 


Y 


J1812 — 1910 









430991 


37 


74 


892 


00 


51804 




55206 


9 


.3 


41 


N 


.T1814-1744 






3 


975905 


744 


70 


792 


00 


51212 




54505 


9 


.0 


50 


P 


J1815— 1738 









198436 


77- 


85 


728 


00 


51157 




55205 


11 


.1 


96 


N 


J1820— 1529 









333243 


37 


91 


772 


00 


51244 




55206 


10 


.8 


41 


N 


J1821 — 1419 






1 


656010 


894 


50 


1123 


00 


51410 




54505 


8 


.5 


53 


N 


J1824— 1945 


B1821- 


19 





189335 


5 


23 


224 


65 


51844 




55206 


9 


.2 


99 


N 


J1825-0935 


B1822- 


09 





769006 


52 


50 


19 


38 


51844 




55073 


8 


.8 


85 


Y 


J1825- 1446 


B1822- 


14 





279187 


22 


68 


357 


00 


51844 




55205 


9 


.2 


96 


N 


J1826- 1334 


B1823- 


13 





101487 


75 


25 


231 


00 


50749 




55530 


13 


.1 


174 


Y 


J1828— 1057 









246328 


20 


70 


245 


00 


51805 




55206 


9 


3 


70 


N 


J1828- 1101 









072052 


14 


81 


607 


40 


51214 




55206 


10 


.9 


42 


N 


J1830-1059 


B1828- 


11 





405043 


60 


03 


161 


50 


51133 




55206 


11 


.1 


212 


N 


J1831-0952 









067267 


8 


32 


247 


00 


51301 




55206 


10 


.7 


77 


N 


J1832-0827 


B1829- 


08 





647293 


63 


88 


300 


87 


51844 




55206 


9 


.2 


80 


N 


J1833-0827 


B1S30- 


08 





085284 


9 


17 


411 


00 


50748 




55206 


12 


.2 


132 


P 


J1834-0731 









512980 


58 


20 


295 


00 


51632 




55206 


9 


.8 


63 


N 


J1835-0643 


B1832- 


06 





305830 


40. 


46 


472 


90 


51529 




55206 


10 


.1 


91 


N 


J1835-1106 









165907 


20 


61 


132 


68 


51945 




55530 


9 


.8 


105 


Y 


J1837-0604 









096294 


45. 


17 


462 


00 


51089 




55206 


11 


.3 


91 


N 


J1838-0549 









235303 


33 


43 


274 


00 


51691 




55206 


9 


.6 


49 


N 


J1839-0905 









418969 


26 


03 


348 


00 


51410 




55206 


10 


.4 


60 


N 


J1841-0524 









445749 


233 


72 


289 


00 


52150 




55507 


9 


.2 


134 


Y 


J1842-0905 









344643 


10 


49 


343 


30 


51460 




55206 


10 


.3 


48 


N 


J1843-0355 









132314 


1 


04 


797 


60 


51159 




55206 


11 


.1 


42 


N 


J1843-0702 









191614 


2 


14 


228 


10 


51692 




55206 


9 


.6 


61 


N 


J1844-0256 









272963 






820 


20 


51559 




55206 


10 


.0 


109 


N 


J1844-0538 


B1841- 


05 





255699 


9 


71 


412 


80 


51844 




55206 


9 


.2 


75 


N 


J1845-0743 









104695 





37 


281 


00 


51633 




55206 


9 


.8 


50 


N 


J1847-0402 


B1844- 


04 





597769 


51 


71 


141 


.98 


51844 




55206 


9 


.2 


75 


N 


J1853-0004 









.101436 




57 


438 


20 


51411 




55206 


10 


.4 


37 


N 


J1853 + 0011 












54 






51148 






11 






N 



1 Y: Glitch detected in this work; N: No glitch detection; P: Previously known glitch before data span. 
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timing solution across the entire data span. In such cases, 
multiple timing solutions were required. 

In order to obtain precise and accurate timing solu- 
tions (including glitch parameters), it is essential to have 
well-determined pulsar positions. For some pulsars, positions 
from the ATNF Pulsar Catalogue were insufficiently accu- 
rate and we therefore determined positions from our data. 
Initially, we obtained timing residuals using the positions 
(and proper motions) provided by the Catalogue. We fitted 
for these parameters using the "Cholesk y" method that ac - 
counts for the effects of correlated noise (I Coles et "all l201l| y 
For each fit, we used the longest data span in which no glitch 
event was observed. The resulting positions were held fixed 
in subsequent processing. 

Glitch events are recognised by a sudden discontinuity 
in the timing residuals relative to a solution based on earlier 
data. For glitches with a fractional size Av^jv larger than 
~ 1CP 6 , the residuals change by a large fraction of the pulse 
period in a few days and phase connection is normally lost 
as illustrated in sub-plot a) of Figure [2] By analysing short 
sections of post-glitch data, phase coherence can normally 
be recovered and an approximate value for the frequency 
glitch Av s determined. As shown in sub-plot b) of Figure 
[2] glitches with smaller fractional sizes, typically ~ 10~ 9 , 
have no loss of phase coherence over several hundred days 
and the post-glitch pulse frequency is easy to determine. 
Even smaller glitches, with Av g /v ~ I0 -10 , are often hard 
to distinguish from irregular timing noise and so the ob- 
served sample of these is incomplete. For each pulsar data 
set exhibiting a possible glitch, we used the glitch plug- 
in of Tempo2 to determine the variations of the pulse fre- 
quency and its first time derivative as a function of time. 
The glitch plug-in realises this by carrying out a sequence 
of local fits for these two parameters to the timing residu- 
als. Typically, we included five or six observations in each fit 
(spanning from about two to six months). After completing 
the sequence of local fits, a list of dates, pulse frequencies 
and pulse frequency derivatives are obtained for each glitch- 
ing pulsar. 

In Tempo2, the additional pulse pha se induced by 
a glit ch is described by Equation (121) in lEdwards et al.l 

g = A<f> + Av p (t - t s ) + ^Av p (t - t s ) 2 + (1) 
[1 - e-C-'^^A^Td 

where the glitch event is modelled by an offset in pulse phase 
A(f> and the permanent increments in the pulse frequency 
Av p and first frequency derivative Av p , in addition to a 
transient frequency increment Av^ which decays exponen- 
tially to zero with a timescale t&. The phase offset Acj> is 
needed to allow for uncertainty in the glitch epoch t s . An 
initial estimate of t s was taken to be halfway between the 
last pre-glitch observation and the first post-glitch observa- 
tion. Initial estimates of Av p and Av p were given by the 
GLITCH plug-in. Improved values were obtained by includ- 
ing the glitch model in the timing model and subsequently 

5 Corrigendum for Equation (121) in.Ejiwards ct al. (200(|): the 
fourth term of the right hand side of the equation should be [1 — 
e-( t r r - t g)A]At/ t r, rather than [1 -e'^-^f^] Ai/ t (^ sr -t s ). 



using Tempo2 to fit for the glitch parameters. For our work, 
we extended the Taylor series in equation {TJ to include AD P 
to characterise the long-term variations in v. These param- 
eters and their corresponding uncertainties were obtained 
from a Tempo2 least-squares-fit to a segment of data typ- 
ically spanning ~200d to ~3000d across the glitch event, 
with the glitch epoch around the centre of the data range. 
The long-term variations in pulse frequency were described 
by a truncated Taylor series, <j>(t) = 4>a + vt + ^vt 2 + \i>1? . 
Fits including Td are more complicated. As Tempo2 imple- 
ments only a linear fitting algorithm, it is necessary to have 
a good initial estimate for Td. The estimate can be realised 
by two steps. In the first step, an estimate for Td was ob- 
tained by eye by inspecting the post-glitch i> variations. In 
the second step, the first-step value was introduced into the 
fitting. By increasing or decreasing Td, one can eventually 
find a Td which minimises the post-fit x 2 . This Td was de- 
termined as the estimate and subsequently included as part 
of the Tempo2 fit. We note that, when a fit included an ex- 
ponential recovery, the post-glitch data range was selected 
to be larger than the recovery timescale, Td. The changes in 
the pulse frequency and its first derivative at the glitch are 
then described as 

Av g = Av p + Av d (2) 

and 

Au s = Au p -^-, (3) 

with their uncertainties obtained using standard error prop- 
agation equations. In addition, a factor Q = Avd/Ai/ g can 
be defined, describing the fraction of glitch recovery. In a 
few cases, after following this procedure the timing residuals 
revealed a shorter-timescale exponential recovery. In these 
cases, the second exponential was fitted, initially holding the 
parameters of the first recovery fixed, and then finally fitting 
for all parameters of both recoveries. 

For some glitches, the glitch epoch could be determined 
by requiring that the pulse phase was continuous over the 
glitch, i.e., that A<f) — 0. However, a unique solution is only 
possible when both the amplitude of the glitch and the in- 
terval between the last pre-glitch observation and the first 
post-glitch observation are small, such that A<f) is less than 
one period between the bounding observations. For situa- 
tions in which this was not possible, we checked the litera- 
ture to determine whether a precise glitch epoch had already 
been published. If so, then we used the published epoch for 
the rest of the analysis. If not, the glitch epoch t s was kept 
at halfway between the last pre-glitch observation and the 
first post-glitch observation, with an uncertainty of half the 
observation gap. To take account of this uncertainty for the 
glitch parameters, we assume a linear dependence on the 
epoch for each of the glitch parameters. The fitting routine 
was carried out again with t g close to the epoch of the first 
post-glitch observation. A difference between the original 
and the new values for each parameter could then be ob- 
tained. The final uncertainty was then the quadrature sum 
of the parameter difference and its original uncertainty. For 
glitches that have a large epoch uncertainty and/or large 
exponential recoveries, the epoch uncertainty term generally 
dominates the final parameter uncertainties. 

Slow glitches are difficult to recognise from timing resid- 
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uals alone and are best identified in plots of v versus time. 
Their identification is somewhat subjective and they can- 
not be fitted with standard glitch analyses. In this paper (in 
i|5.3[) we describe slow glitches detected in one pulsar, PSR 
J1539-5626. 



4 RESULTS 

The data sets for the 165 pulsars in the sample were pro- 
cessed, and 36 pulsars were observed to have glitched (indi- 
cated with a "Y" in the last column in Table Q]). A total of 
107 glitches were detected, among which 46 are new detec- 
tions. We identified exponential recoveries for 27 glitches. 
A total of 22 previously published glitches are within our 
data span, but we were unable to identify these events. This 
is mainly because the sampling of our observations is often 
insufficient, such that glitches with a small fractional size 
(Av s /v < 10~ 9 ) are hard to detect. For the same reason, 
only those exponential recoveries with a timescale between 
a few tens to a few hundred days are detectable; any expo- 
nential recoveries with a timescale shorter than a few tens 
of days are likely to have been missed. 

Table[2]gives the positions and proper motions in J2000 
coordinate for each glitching pulsar. The positions for 28 
pulsars are from the ATNF Pulsar Catalogue. As described 
in §3, we fit for the positions for a further eight pulsars. All 
of the proper motions are from the ATNF Pulsar Catalogue. 

Table IA1I lists the pre-, inter- and post-glitch timing 
solutions for the glitching pulsars. For each pulsar, the ta- 
ble contains the pulsar name, the interval relative to glitch 
number, v, i>, i>, reference epoch, fitted data span, number of 
ToAs, post-fit rms residuals, reduced \ 2 an< A the number of 
degrees of freedom for the least-squares fit. For pulsars with 
exponential post-glitch recoveries, to avoid contaminating 
the long-term post-glitch parameters, the start of post-glitch 
data span is at least two decay timescales from the glitch. 
These solutions include long-term timing noise and so are 
only valid within the fitted data range; they cannot be used 
for extrapolation. 

Table IA2I contains the parameters for each observed 
glitch. The second column gives a reference number for each 
glitch and the glitch epochs are given in the third column. 
The fourth column indicates whether the glitch is new (N) 
or has been previously published (P). References for previ- 
ously published glitches may be found in the web databases. 
For each glitch parameter, we give two uncertainties. The 
Tempo2 1 a uncertainties are given in the first pair of paren- 
theses. If inclusion of effect of the glitch epoch uncertainty 
made a significant difference, the final uncertainty is given 
in the second pair of parentheses. Note that errors refer to 
the last digit quoted. The number of observations, the fitted 
data span, the post-fit rms residuals and the reduced \ 2 an d 
degrees of freedom are listed in columns 11, 12, 13 and 14, 
respectively. 

In each of the sub-sections below, we describe the ob- 
served glitch events for each pulsar in more detail. In Figures 
[31 to 1111 for the 36 glitching pulsars, we show the evolution 
of pulse frequency and its first time derivative within our 
data span. For convenience, observed glitches are numbered 
as in Table El 



4.1 PSR J0729 1448 

A data gap lasting for ~ 6yr exists in the data set of 
this pulsar. During our data span of the recent three 
years, this pulsar exhibited four glitches (see Figure [3]). 
The first three glitches were small (Av E /v ~ 10 -8 ). The 
fourth glitch was significantly larger {Av R jv ~ 6 x 10~ 6 ). 
These glitch events have been repor ted bv lWeltevrede et al.l 
|20ld) and lEspinoza et al.l (|2011al l. The Parkes data are 
unfortunately not well sampled. The thre e small glitches 
were id entified with prior knowledge from lEspinoza et al.l 
|2011al ). For the same reason, it is impossible to evaluate 
the permanent change in ii and other long-term parame- 
ters. IWeltevrede et all ()2010h reported the large glitch that 
occurred at MJD 54711(21). Our analysis provides a more 
precise epoch o f MJD 54681(9) which is consistent with 
lEspinoza etafl (|2011al ) result of MJD 54687(3). 



4.2 PSR J0742 2822 (PSR B0740-28) 

In tota l, seven glitch events have been reported for this 
pulsar (jD'Alessandro et al.l ll993; Janss en fc Stappersll2006l ; 
lEspinoza et al.l 201 lal ). In Figure |3j we present our 17-yr 
data span. No new glitches were detected. A glitch at MJD 
~ 55020 can clearly be seen. However, we were unable to 
detect the four small previously reported glitches covered 
by the data set. For the observed glitch, Es pinoza et al.l 
(|2011aT ) gave Au/u — —0.372(96), corresponding to Av — 
225(58) x 10 -15 s -2 . Our measurement of v presented in Fig- 
ureEJand Table IMl shows Av = -1.3(3) x 10" 15 s" 2 , despite 
the evident noise. 



4.3 PSR J0834 4159 

This pulsar was not previously known to glitch, but we iden- 
tify a small glitch at MJD ~ 53415. Figure [3] shows the 10-yr 
evolution of pulse frequency and pulse-frequency derivative 
of this source observed at Parkes. Both of the measured Av 
and v exhibit noise, v shows a small permanent change at 
the glitch event (Table |A"2)) . Our observations do not reveal 
any post-glitch relaxation process. 



4.4 PSR J0835 4510 (PSR B0833-45) 

The Vela pulsar has undergone 16 known glitch events over 
a period of ~ 38 yr (for a complete list of these glitches, 
see the ATNF Pulsar Catalogue glitch table or the Jodrell 
Bank Glitch Catalogue). Thirteen have a fractional glitch 
size larger than 10 -6 . In Figure [3] we present the varia- 
tions of pulse frequency and its first derivative spanning the 
last ~ 15 yr. Four glitches wer e detected . Thes e events have 
been r e ported and anal y sed by Flanagan (1996); Wan g~et al.l 
||2000D; iDe Luca et all (11999); iDodson et all l|2002l . I2004D 
and lFlanagan fc Buchnerl ( 2006T ). 

As Figure [3] shows, these glitches are large, with 
Avgjv > 2 x 10~ 6 . Each of the post-glitch behaviours ex- 
hibits both exponential and linear recoveries. We attempted 
to model each of the glitches including both of the exponen- 
tial and linear recoveries. Each of the glitches is discussed in 
more detail as b elow. 

For glitch 1, 1 Wang et aD (|2000h reported an exponential 
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Table 2. Position and proper motion parameters for 36 glitching pulsars. 



PSR J R. A. Dec. Position epoch fi a fig References 
(h:m:s) (° ' ") (MJD) (mas yr" 1 ) (mas yr" 1 ) 



J0729- 


-1448 




— 14-48-^fi 


51367 






1 


J0742- 


2822 


07:42:49.058(2) 


-28:22:43.76(4) 


49326 


-29(2) 


4(2) 


2,3 


J0834- 


-4159 


08:34:17.815(8) 


-41:59:36.01(9) 


52347 


- 


- 


This work 


J0835- 


-4510 


08:35:20.61149(2) 


-45:10:34.8751(3) 


51544 


-49.68(6) 


29.9(1) 


4 


J0905- 


5127 


09:05:51.94(5) 


-51:27:54.0(4) 


54072 


- 


- 


This work 


J1016- 


5857 


10-16-21 16fPl 


— 58-57-1 2 1 ("1 "1 


52717 






5 


J1048- 


-5832 


10:48:12.2(1) 


-58:32:05.8(8) 


50889 


- 


- 


6 


J1052- 


5954 


10:52:38.11(7) 


-59:54:44.1(5) 


51683 


- 


- 


7 


J1105- 


6107 


11:05:26.17(4) 


-61:07:51.4(3) 


50794 


- 


- 


6 


J1112- 


6103 


11:12:14.81(4) 


-61:03:31.1(6) 


51055 


- 


- 


8 


J1119- 


6127 


1 1 -19-14 30f2^ 


— 61 -27-49 5l'2') 


51485 






9 


J1301- 


6305 


13:01:45.76(14) 


-63:05:33.9(12) 


51206 


- 


- 


8 


J1341- 


6220 


13:41:42.63(8) 


-62:20:20.7(5) 


50859 


- 


- 


6 


J1412- 


6145 


14:12:07.69(5) 


-61:45:28.8(6) 


51186 


- 


- 


8 


J1413- 


6141 


14:13:09.87(9) 


-61:41:13(1) 


51500 


- 


- 


7 


J1420- 


6048 


14-20-08 237(1 6) 


— 60-48-16 43C15'l 


51600 






10 


J1452- 


6036 


14:52:51.898(8) 


-60:36:31.35(6) 


51630 


- 


- 


7 


J1453- 


6413 


14:53:32.684(8) 


-64:13:15.81(7) 


52608 


-16(1) 


-21.3(8) 


This work 


J1531- 


5610 


15:31:27.91(1) 


-56:10:55.0(1) 


51448 


- 


- 


7 


J1614- 


5048 


16:14:11.29(3) 


-50:48:03.5(5) 


50853 


- 


- 


6 


J1646- 


4346 




—4^-4^-48^ 


52792 






r 1 nio ~\\TC\v\i~ 
_L HIS WUlA 


J1702- 


4310 


17:02:26.94(5) 


-43:10:40(2) 


51597 


- 


- 


7 


J1709- 


4429 


17:09:42.728(2) 


-44:29:08.24(6) 


50042 






6 


J1718- 


3825 


17:18:13.565(4) 


-38:25:18.06(15) 


51184 






8 


J1730- 


3350 


17:30:32.28(6) 


-33:50:28(4) 


53826 






This work 


J1731- 


4744 


17:31:42.17(7) 


-47:44:37(2) 


54548 






This work 


J1737- 


3137 


17:37:04.29(4) 


-31:37:21(3) 


51234 






1 


J1740- 


-3015 


17:40:33.82(1) 


-30:15:43.5(2) 


52200 






3 


J1801- 


2304 


18:01:19.829(9) 


-23:04:44.2(2) 


50809 






11 


J1801- 


2451 


18:01:00.016(8) 


-24:51:27.5(2) 


53348 


-11(9) 


-1(15) 


12 


J1803- 


2137 


18:03:51.4105(10) 


-21:37:07.351(10) 


51544 


11.6(18) 


14.8(23) 


13 


J1809- 


1917 


18:09:43.132(6) 


-19:17:40(1) 


54632 






This work 


J1825- 


0935 


18:25:30.629(6) 


-09:35:22.3(3) 


53300 


-13(11) 


-9(5) 


14,3 


J1826- 


-1334 


18:26:13.175(3) 


-13:34:46.8(1) 


52400 


23.0(25) 


-3.9(31) 


14,15 


J1835- 


1106 


18:35:18.41(7) 


-11:06:15(4) 


53882 






This work 


J1841- 


0524 


18:41:49.32(5) 


-05:24:29.5(12) 


52360 






5 



er m otions: 1 -iMorris et alj <2002h : 2 
5 - iHobbs et alj j2004h :6 



Refere nces f or positions and pro 
Jl997h; 4 -iDodson et al.l ||2003| _ 

- [Manchester et al l ll 200lh; 9 - ICamilo et al.l feOOOh; 10 - iD'Amico et alj (l200lh: 11 - iFrail et al.l jl993h : 12 
IZeiger et al.l j2008h; 13 -iBrisken et al.l j20odT: 14 -lYuan et alj fcOlOal); 15 -iPavlov et alj <2od3i 



IHobbs et al.l II2004I): 3 -iFomalont et al.l 
Wang et a.1.1 d200Ctl : 7 - iKramer et alj ll2003h : 8 



recovery with a time constant 916(48) d. However, the post- 
event v variations shown in Figure indicate that the expo- 
nential recovery completes within ~ 200 d and the long-term 
evolution exhibits a linear recovery. Fitting the timing phase 
residuals showed that the exponential timescale is 186(12) d, 
with Q = 0.030(4) (Tabl e lAH). Glitch 2 was captured with 
high time resolution by Dods on et all l|2002i ). Four short- 
term exponential decays were identified with the smallest 
timescale just ~ 1.2 min. Parkes data, however, are not suf- 
ficient to resolve these short-term recoveries. However an ex- 



ponential recovery that completes in ~ 100 d can be seen in 
our data. Fitting gave Q ~ 0.02 and ~ 125 d. For glitch 
3, the exponential recovery is characterised by Q ~ 0.009 
and Td ~ 37 d. The most recent glitch 4 also exhibits an 
exponential recovery. Our fitting showed Q ~ 0.0119 and 
r d ~ 73 d. 

At least in the long term, the post-glitch behaviour is 
dominated by the linear recovery of v. This is superimposed 
on the shorter-term exponential decays and persists until 
the next glitch. As Table [ATI shows, the observed values of 
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Figure 3. Glitches in PSRs J0729-1448, J0742-2822, J0834-4159 and J0835-4510. The top panel shows the pulse-frequency residuals 
Aw, obtained by subtracting the (extrapolated) pulse frequency derived from the v and i> values of the first pre-glitch solution; the middle 
panel is an expanded plot of Av where the mean of each inter-glitch (or post-glitch) solution is subtracted; the bottom panel shows the 
variations of the pulse-frequency first time derivative v. The glitch epochs are indicated by vertical dashed lines, and the numbers at the 
top of each dashed line denote the sequence of glitches detected within our data span. 



v representing the slope of this long-term linear recovery- 
are relatively large for the Vela pulsar and also that they 
change signficantly after each glitch. Fitting for Ai> p along 
with the other glitch parameters was generally difficult. The 
superimposed timing noise and the shorter data spans used 
for the glitch fitting often led to values somewhat different to 
those obtained by differencing the long-term fits for v given 
m Table ED or to an insignficant value. For example, glitch 
3 has a fitted value of 304(23) x 10" 24 s -3 but the difference 
between the post- and pre-glitch values of v in Table I All is 
491(7) x 10" 24 s -3 . 



4.5 PSR J0905 5127 

No glitch events have previously been reported for this pul- 
sar. Figure [4] presents the evolution of v and v observed at 
Parkes. The entire data span is ~ 16 yr, but there exists 
a data gap lasting for ~ 4yr. Two glitch events were de- 
tected. Both are small, with a fractional glitch size ~ 10~ 8 . 
The available observations are not sufficient to study the 
post-glitch behaviour for glitch 1. For glitch 2, no signifi- 
cant post-glitch recovery was observed. 



4.6 PSR J1016 5857 

In Figure 2] the variations of v and i> of this pulsar for 
~ 11 yr are shown. Two glitches were detected. They are 
similar with Av g /i> ~ 2 x 10~ 6 and Az> g //> ~ 4 x 10~ 3 . 
The different slopes of v before and after glitch 1 imply a 
permanent change in v\ fitting showed that Ai> p = 69(7) x 
10 -24 s -3 , approximately consistent with the v values in 
Table [ATI For glitch 2, the available data are not sufficient 
to characterise the long-term post-glitch relaxations. 



4.7 PSR J1048 5832 (PSR B1046 -58) 

For PSR J1048— 5832, the evolu tion of v a n d v span- 
ning 20 yr is shown in Figure H IWang et all (|2000l ) and 
lUramal (|2002|) have publis hed details for glitches 1, 2 and 
3. IWeltevrede et all |20ld ) discovered glitch 6. We report 
here glitches 4 and 5 as new discoveries. Glitch 4 is large 
with Av s /v ~ 1.8 x 10~ 6 , whereas glitch 5 is much smaller 
with a fractional size ~ 2.5 x 10~ 8 . As shown in Figure [4l 
for these two glitches there is little evidence for exponential 

recoveries. 

For glitches 2 and 3. IWang et all (|2000h included expo- 
nential terms to model the post-glitch behaviour; the time 
constants were assumed to be 100 d and 400 d, respectively. 
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Figure 4. Glitches in PSRs J0905-5127, J1016-5857, J1048-5832 and J1052-5954. See Figure|3]for a description of each sub-plot. 



lUramal l|2002f ) observed glitch 3 with high observing cadence. 
Two exponential recoveries were detected; the timescales are 
32(9) d and 130(40) d, respectively. As shown in Figure [4] 
for both of the glitches 2 and 3, the post-glitch v variations 
exhibit significant noise. For glitch 2, our fitting for the ex- 
ponential recovery showed Q = 0.026(6) and r d = 160(43) d. 
We note that, for this glitch, because there is only one pre- 
glitch measurement of z>, so we are unable to measure the 
permanent change in v. For glitch 3, our fitting showed 
Q = 0.008(3) and r d = 60(20) d. Table IAT1 shows signficant 
values of v for all except glitch 2, with significant variations 
from glitch to glitch. On average, the values are about an 
order of magnitude smaller than those for the Vela pulsar. 

4.8 PSR J1052-5954 

The available data set for this pulsar contains a data gap 
of ~ 6yr. Figure |3] shows the evolution of v and v after the 
da ta gap. The detect e d glit ch at MJD ~ 54495 was reported 
by I Weltevrede et alj (|2010T ). An exponential relaxation and 
a significant permanent increase in spin-down rate \v\ can 
be identified in the post-glitch data. Fitting to the timing 
residuals indicated that ~ 0.067 of the glitch recovered in 
~ 46 d. 

4.9 PSR J1105 6107 

Three g litch events have p reviously been identifie d for this 
source (|Wang et al.l |2000| ; IWeltevrede et alj|20ld ). In Fig- 



ure [5] we present the evolution of v and v for ~ 16 yr of this 
pu lsar. We confirm the previously detected g l itches 1 and 
3 dWang et all |2000| ; IWeltevrede et all l2010h . IWang et all 
l|2000! ) reported a small glitch occurred at MJD - 50610. 
However, we found that the timing behaviour of this source 
around this epoch is more likely to be dominated by timing 
noise. We report new glitch events as glitches 2 and 4. As 
shown in Figure [5j the post-glitch behaviour is noisy and no 
exponential recoveries were observed. There appears to be a 
persistent increase in \v\ at the time of each of glitches 1, 2 
and 3. For glitch 4, the available data span is not adequate 
to study the post-glitch behaviour. 



4.10 PSR J1112 6103 

As shown in Figure [5] two large glitch events were detected 
in this pulsar, with Av g /u ~ 10~ 6 . For glitch 1, the ob- 
served variations of v indicate a large change in i>; fitting 
gave Ai/p ~ 240 x 10~ 24 s -3 . No exponential recovery was 
observed for this glitch. For glitch 2, a long-term exponential 
relaxation was observed, which is characterised by Q ~ 0.022 
and rd ~ 300 d. 



4.11 PSR J1119 6127 

For PSR J1119-6127, Figure |5] shows the evolution of the 
observed pulse frequency and its first derivative spanning 
~ 13 yr. Three glitches were observed. The first is small 
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Figure 5. Glitches in PSRs J1105-6107, J1112-6103, J1119-6127 and J1301-6305. See Figure|3]for a description of each sub-plot. 



with AUfr/u ~ 4 x 10 9 as was reported by ICamilo et al,l 
(2000). The second and t hird glitches are much l arger and 
were studied in detail bv IWeltevrede et al.l l|201ll 'l. Our re- 
sults are generally consistent with theirs. 



4.12 PSR J1301 6305 

Figure [5] shows the evolution of v and v for PSR J1301— 6305 
over ~ 11 yr. We detected two large glitch events. Glitch 
1 has Afg/f ~ 4.6 x 10 -6 and for glitch 2 the fractional 
frequency change is about half this. For glitch 1, an expo- 
nential recovery was identified; fitting gave Q ~ 0.0049 and 
Td ~ 58 d. As shown in Figure [5] the pre- and post-glitch in- 
tervals show clear linear recoveries. As Table IAT1 indicates, 
the i> values are all ~ 250 x 10~ 24 s~ 3 . Because of tim- 
ing noise, we were unable to fit for the Ai> p changes at the 
glitches. 

4.13 PSR J1341 6220 (PSR B1338 62) 

This pulsar is we ll kn own to have freq u ent g litches - 
I Wang et all (|2000h and IWeltevrede et all (|2010h have re- 
ported 14 glitches. Figure [5] shows the evolution of v and v 
for ~ 16 yr, where a total of 17 glitches are presented. We re- 
port here the new detections of nine glitch events. For glitch 
6, an exponential decay with Q ~ 0.0112 and Td ~ 24 d was 
detected. Unfortunately, for the other glitches the observa- 
tions are insufficient to study the post-glitch behaviour. 



4.14 PSR J1412-6145 

PSR J1412— 6145 has not previously been known to glitch. 
Here, we report the discovery of a large glitch with A^ g /^ ~ 
7.2 x 10~ 6 that occurred at MJD ~ 51868. As shown in Fig- 
ure |6j there was a clear exponential recovery with timescale 
~ 60 d, a significant increase in |z>| at the time of the glitch 
and a slow linear recovery of part of this increase. 

4.15 PSR J1413-6141 

Figure [6] presents seven new glitch events detected in this 
pulsar over a 12.6-yr data span. Among these events, three 
are small (Ai> s /v ~ 10 -8 ), while the other four are larger, 
with a fractional size > 10 -6 . Exponential post-glitch recov- 
eries are not observed for these glitches. Significant values 
of v are seen after each of the latest four glitches (cf. Ta- 
ble lAip . We were able to fit for Ai> p for glitch 4, giving 
a value of 491(42) x 10~ 24 s~ 3 ; this is consistent with the 
difference between the post- and pre-glitch solutions for v, 
which is 457(40) x 10" 24 s -3 . 

4.16 PSR J1420 6048 

Figure [6] shows that five glitch events were observed in 
this p ulsar. Glitch 4 was first reported bv IWeltevrede et al.l 
(2010). All of these glitches are large, with Av g /v ~ 10" 6 . 
The post-glitch slow-down rates exhibit linear decays, and 
changes in v are observed (Tables [All and IA2|) . For glitch 
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Figure 6. Glitches in PSRs J1341-6220, J1412-6145, J1413-6141 and J1420-6048. See Figure|3]for a description of each sub-plot. 



2, an exponential recovery was measured, with obtaining 
Q ~ 0.008 and Td ~ 99 d. For glitch 5, the available data are 
not sufficient to study the post-glitch behaviour. 



4.17 PSR J1452 6036 



exponential recovery, an offset in v at the time of the glitch 
and a slow linear recovery. Our fitting of the exponential 
term showed that ~ 0.007 of the glitch recovered within a 



timescale of < 
(Table El"]) . 



76 d and the long-term v is ~ 20 x 10 s 



This pulsar was not previously know n to glitch. Since th e 
end of the Multibeam Survey timing l|Kramer et al.l f2003), 
no observations were made until the start of the Fermi 
project. Hence, a data gap lasting for ~ 5yr exists. Fig- 
ure [7] shows the evolution of v and v. A small glitch event 
with Az/g/z/ ~3x 10~ 8 was detected at MJD ~ 55055. The 
available data are not adequate to study the post-glitch be- 
haviour. 

4.18 PSR J1453-6413 (PSR B1449-64) 

No glitch event has previously been reported for this pulsar. 
Figure [7] shows the evolution of v and v for ~ 12 yr. We 
detected a very small glitch with Au g /u ~ 3 x 10 -10 . We 
cannot comment on the post-glitch behaviour since the data 
are insufficient. 



4.20 PSR J1614 5048 (PSR B1610 50) 

PSR J1614-5048 has been observed at Parkes for ~ 20 yr. 
As shown in Figure [7] two glitches were detected. Both of 
the events are large, with At/ g /u > 6 x 10~ 6 . This pulsar ex- 
hibits remarkable timing noise; the large-scale fluctuations 
in i> reflect this. As a result, phase-connected timing resid- 
uals cannot be obtained for the entire data range between 
the two glitch events. We thus report the timing solutions 
for this data span in two sectio ns (see Tabl e IA1I) . Glitch 1 
has previously been reported by I Wang etail (2000); our re- 
sults for Av s /v and Az> g /z> are consistent with theirs. Glitch 
2 is a new detection. Values of v given in Table rM] show 
significant variations, but these are likely to be contami- 
nated by the timing noise. Despite the noise, there does ap- 
pear to be a signficant linear recovery after glitch 1 with 
v ~ 200 x 10~ 24 s' 3 . 



4.19 PSR J1531 5610 

PSR J1531 — 5610 was not previously known to glitch. Figure 
[7] shows a large glitch event at MJD ~ 51730, detected by 
Parkes timing. As in PSR J1412— 6145, this glitch has an 



4.21 PSR J1646 4346 (PSR B1643 43) 

Figure [8] presents the evolution of v and v for this pulsar 
for ~ 16 yr. A glitch event was detected at MJD ~ 53875. 
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Figure 7. Glitches in PSRs J1452-6036, J1453-6413, J1531-5610 and J1614-5048. See Figure|3]for a description of each sub-plot. 



This is the first reported glitch for this pulsar and it has a 
fractional size ~ 8.8 x 10~ 7 . There is a clear linear recovery 
from a presumed earlier glitch before the observed glitch. 



4.22 PSR J1702 4310 

This pulsar was not previously known to glitch. Here, as 
shown in Figure [8] we report our discovery of a glitch event. 
Our observations suggest both exponential and linear decays 
following the glitch and a linear decay with almost the same 
slope preceding the glitch. Fitting for the exponential decay 
gave Q ~ 0.023 and Td ~ 96 d. A data gap of ~ 3yr exists 
in the pre-glitch data span. However, there is little period 
noise and no phase ambiguity across this gap. 



4.24 PSR J1718 3825 

This pulsar was not previously reported to glitch. Figure 
[H] shows the evolution of v and v for ^13 yr. A glitch was 
detected at MJD ~ 54910. This event is small with Au s /u ~ 
2 x 10~ 9 . 



4.25 PSR J1730 3350 (PSR B1727 33) 

Two glitch events that occurred at MJDs ~ 48000 
and ~ 52107 were p reviously detected for this pulsar 
l|johnston et all 1 19951 : lEspinoza et al.1 ^Ollal ). Both are 
large, with Avjv > 3 x 10" 6 . Figure [9] shows the MJD 
~ 52017 glitch. There is a clear linear recovery following the 
glitch and a small exponential recovery with Q ~ 0.01 and 
timescale of ~ 100 d. 



4.23 PSR J1709-4429 (PSR B1706 44) 

In Figure [8] we present the evolution of v and v of this pul- 
sar for ~ 20 yr. Four glitches were d etected. These glitches 
are large, with Au R /v > 1 x 10" 6 . iJohnston et alj (|l995h 
and IWeltevrede et al.l (|2010l ) have reported glitches 1 and 
4, but no post-glitch recoveries were reported. As Figure [8] 
shows, all four glitches show significant post-glitch recov- 
eries, most with both exponential and linear components. 
Dramatic slope changes in the linear recoveries are seen af- 
ter each glitch (Table lAlfl . Just a small fraction of each 
glitch recovers exponentially, with time constants ~ 100 d 
(Table O. 



4.26 PSR J 1731 4744 (PSR B1727 47) 

Three glitches have been repo rte d for this pulsar 
ilD'Alessandro fc McCullochl Il997l ; IWang et al.l l200d : 
lEspinoza et al.l 201 lal ). Figure [9] shows these numbered 1, 
2 and 3 and a newly detected glitch 4. Glitches 1 and 3 
have Au g /u ~ 1 X 10 -7 , whereas glitches 2 and 4 have 
Av g /v ~3x 10~ 9 . For glitch 3, we measured an exponential 
recovery, with obtaining Q ~ 0.073 and Td ~ 210 d. There 
is so me evidence for lin ear recoveries following each glitch 
(cf., IWang et al.1 l2000h . The measured values of v are 
small, ranging from 1 x 10" 24 s" 3 to 20 x 10" 24 s" 3 . For 
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Figure 8. Glitches in PSRs J1646-4346, J1702-4310, J1709-4429 and J1718-3825. See Figure|3]for a description of each sub-plot. 



glitch 2, IWang et al.l l|2000l ) reported the glitch epoch as 
MJD 50703(5). We obtained a more accurate epoch MJD 
50715.9(8) by solving for phase continuity across the glitch. 



4.27 PSR J1737 3137 

Since the end of the Multibeam Pulsar Survey timing 
l|Morris et all 12002? ) . this source was not observed until the 
commencement of the Fermi project, leaving a data gap of 
~ 7yr. Three glitches have been reported f or this pulsar 
l|Weltevrede et al.ll2010l ; lEspinoza et alj2011al ). We detected 
the most recent glitch (Figure |SJ. This glitch is large, with 
Avg/v ~ 10 -6 . There is evidence for a signficant change in v 
at the time of the glitch, and maybe exponential and linear 
recoveries, but the available data are insufficient to be sure. 



4.28 PSR J1740 3015 (PSR B1737 30) 

PSR J1740— 3015 is one of the most frequently glitching pul- 
sars known. Previously, a total of 31 glitches were detected 
over a 25-yr data span. Figure [9] presents the evolution of 
v and v for ~ 13 yr from Parkes timing. During this pe- 



(lUramal 12002; 


Krawczvk et al.l 120031; Ijanssen & Staoperc 


2006: 


Zou et al 


2008 


; Yuan et al. 2010al; Weltevrede et al. 


201C; 


Esphioza et al. 


|2011a[). However, because of the fre- 



four have been published. Glitch 5 is a new detection; it has 
a fractional size of ~ 2.7 x 10~ 6 , which is the largest ever 
seen in this pulsar. Linear recoveries are seen after most of 
the large glitches, including glitch 5, but we have no evidence 
for exponential recoveries. If such recoveries exist, they must 
have timescales of less than a few tens of days. 



4.29 PSR J1801-2304 (PSR B1758 23) 

PSR J1801— 2304 has been f ound to have suffered nine 
glitches in ~ 24 y r dKaspi et al 



1119931; IShemar fc Lvnelll996 
I Wang et afll2000l; Krawczyk et al.l 120031; lYuan et al.ll2ofoa 



quent glitching and the relatively poor sampling of our ob- 
servations, only five major glitches can be detected. The first 



Espinoza et al.ll2011al ). Figure [ToT presents the variations of 
v and v for ~ 20 yr, where eight glitches are shown, the 
last being a new discovery. This is a s mall glitch, with a 
fractional size ~4x 10~ 9 . For glitch 4, IWang et al] (|200dh 
fit it with an exponential recovery, assuming the timescale 
to be 100 d. However, as shown in Figure 1101 our observa- 
tions suggest that a linear recovery with v ~ 40 x 10 - 24 s -3 
(Table lAlj) dominates the post-glitch behaviour. 



4.30 PSR J1801 2451 (PSR B1757-24) 

Figure [10] shows the evolution of v and v of this pulsar 
for ~ 18 yr, where five glitch es are shown. These glitches 
have been reported previously dLvne et al 1 19961; I Wang et all 
2000| IWeltevrede et aUl20ld : lEspinoza et~aHl2011a| ) . Figure 
[10] suggests that the post-glitch recoveries are dominated by 
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Figure 9. Glitches in PSRs J1730-3350, J1731-4744, J1737-3137 and J1740-3015. See Figure|3]for a description of each sub-plot. 



linea r recoveries but there may be expon ential recoveries as 
well |Lvne et al.lll996l ; IWang et all 120001 ) . In this work, we 
fit exponential recoveries for glitches 3 and 5. For glitch 3, 
we obtained Q ~ 0.025 and ~ 200 d and for glitch 5, 
Q ~ 0.0065 and r d ~ 25 d. 

4.31 PSR J1803-2137 (PSR B1800-21) 

PSR J1803-2137 is very similar to PSR J1801-2451 in its 
timing properties with a characteristic age of 15 kyr. Four 
glitches have been detected in this pulsar llShemar fc Lvne] 
19961; I Wang et al.l [200ol; iKrawczvk et al.ll2003l ; lYuan etai] 
2010al ; lEspinoza et all ^Ollal ). Figure [10] shows the two 
glitches detected in this work. Both of these show clear long- 
term linear recoveries with a very similar slope (Tabic lAll) . 
Clear exponential recoveries with sh orter timescales are also 
seen. For glitch 1. 1 Wang et all (|200Ch fitted a short timescale 
(~ 18 d) recovery and a ~ 850 d one to the longer-term de- 
cay. We found evidence in the phase residuals for two short- 
term decays with timescales of ~ 12 d and ~ 69 d; further 
rel axation is dominat ed by the linear recovery. For glitch 
2, lYuan et al.l (|2010aT ) fit an exponential decay, obtaining 
Q = 0.009(2) and r d = 120(20) d. We obtained consistent 
results as Q = 0.00630(16) and r d = 133(11) d. 

4.32 PSR J1809-1917 

This pulsar has previ ously been reported to glitch once 
|Espinoza et al.l l2011al ) . Figure \W\ shows the evolution of 



v and v for ~ 12yr. The glitch is large, with Av s /u ~ 1.6 x 
10~ 6 . Both exponential and linear recoveries can clearly be 
seen in our observations. Fitting showed that the exponen- 
tial recovery is characterised by Q ~ 0.006 and r d ~ 125 d. 
A linear recovery, presumably from an earlier glitch, is also 
seen before the glitch. The recovery after the observed glitch 
is slightly less steep, but both have v ~ 35 x 10~ 24 s -3 . 



4.33 PSR J1825 0935 (PSR B1822 09) 

So far, eight glitches have been detect ed for this pul- 
sar; six were identified as "slow" gl itches JZou et al. 112004 



Shabanova 



12007 



Yuan et alj|2010al ; lEspinoza et alj|2011a 



Lvne et al.l (|2010r ) have shown that these "slow" glitches 
probably are a manifestation of "two-state" magnetospheric 
switching. Figure [TT] presents our observations of the evolu- 
tion of v and v for ~ 9yr. During this pe riod, this pulsar 
was found to undergo three slow g litches |Zou et al.ll2004 
IShabanovall2007l ; lYuan et alj|2010al ). Because of insufficient 
observations, these slow glitches were missed in our data. 
By using the available data, we were able to identify a "nor- 
mal" glitch with Ai> g /i> ~ 1.3 x 10~ 7 . However, a data gap of 
~ 500 d exists shortly after this glitch, making it impossible 
to study the post-glitch behaviour. Th is glitch was observ ed 
in more detail by the Xinjiang group (|Yuan et al.ll2010ah . 
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Figure 10. Glitches in PSRs J1801-2304, J1801-2451, J1803-2137 and J1809-1917. See Figure|3]for a description of each sub-plot. 



4.34 PSR J1826 1334 (PSR B1823 13) 

Five glitch events have been reported for this pulsar 
llShemar fc Lvnd Il996l ; lYuan et all l2010al ; lEspinoza et~ai1 
l2011al ). Figure [TTJ presents the evolution of v and ii for 
~ 13 yr. For glitch 1, our last pre-glitch observation was 
at MJD ~ 53186, and the first post-glit ch observation was 
at MJD ~ 53279. lEspinoza et alj (|2011bJ ) showed that there 
are actually two small glitches in this interval. One was at 
MJD ~ 53206, with Avjv ~ 0.6 x 10~ 9 and the other was 
at MJD ~ 53259, with Av g /v ~ 3 x 10~ 9 . The much l arger 
glitch 2 was also obser ved by both lYuan et all (|2010ah and 
lEspinoza et al.l l|2011aT ). We detected an exponential recov- 
ery for this glitch, characterised by Q ~ 0.007 and r d ~ 80 d. 



4.35 PSR J1835-1106 

Figure [TT] shows the evolution of v and v for 
PSR J1835-1106 over about 10 yr. A small glitch with 
Av z jv ~ 1.6 x 10~" 8 was detect ed at MJD ~ 5 2220. 
This event wa s also observed by IZou et ail |2004 ) and 
lEspinoza et al.l (|2011al ) . As shown in Figure 111! the post- 
glitch frequency residuals exhibit a cubic structure, indi- 
cating a measurable v; fitting showed that this term is 
1.58(13) x 10~ 31 s -4 . These signficant higher-order frequency 
derivatives indicate the presence of noise processes in the 
pulsar rotation. It is possible that this pulsar also has a 
tw o-state magnetos pheric modulation affecting the value of 
v jLvne et alj|20ld ). 



4.36 PSR J1841 0524 

This pulsar ha s been found to undergo three glitches 
ijEspinoza et al.l l2011al 'l. The first two are small, with 
^■ v zl v ^ 2 x 10~ 9 . The latest one is large, with Av g /v ~ 
10~ 6 . Parkes data have a gap, spanning from MJD ~ 52570 
to MJD ~ 53619. As a result, the MJD ~ 53562 glitch was 
missed. Figure [TTJ presents the evolution of v and ii for the 
past five years. At least for glitch 2, the post-glitch behaviour 
exhibits a linear recovery. 



5 DISCUSSION 

We have searched for glitch events in the timing residuals of 
165 pulsars. Out of these, 107 glitches were identified with 
46 new discoveries. Most of these new discoveries occur for 
southern-hemisphere pulsars that cannot be observed by the 
long-term monitoring programmes carried out in the north- 
ern hemisphere. Because of the relatively large gaps between 
many of our observations, there are 22 known glitches that 
are undetectable in our data sets. These missed glitches gen- 
erally have fractional sizes Af g /f between 10 -10 and 10~ 9 . 
The measurement of Av g is also very dependent on the data 
sampling as short-term transients can easily be missed or un- 
derestimated and in some cases our results differ from those 
in the literature. 

In general, the post-glitch behaviour shows two types of 
recovery: a short-term exponential recovery (characterised 
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Figure 11. Glitches in PSRs J1825-0935, J1826-1334, J1835-1106 and J1841-0524. See Figure|3]for a description of each sub-plot. 



by Q and Td) and a longer-term linear recovery in v (char- 
acterised by D). Both can be identified by inspecting the evo- 
lution of v. In most cases, we have insufficient observations 
to study any exponential recoveries with timescale <20 d. 
Most such short-term exponential recoveries will have been 
missed. However, we found that 27 glitches do show mea- 
surable post-glitch exponential recovery with time constants 
between 12 d and 300 d. For more than 90% of the observed 
glitches, values of post-glitch i>, indicating a linear recovery 
in 0, were obtained. For many of these, the v slope was sim- 
ilar before and after the glitch, so the value of AD P was not 
significant. For 13 glitches, the slope change was larger and 
a significant value for Ai> p was obtained. 

The discovery of 46 new glitches allows further study on 
the distribution of glitches and their post-glitch behaviour. 
The discussion on these two aspects is presented in i]5.1l and 
i )5.2l respectively. 



5.1 The distribution of glitches 

In Figure 1121 the upper panel is a histogram of the frac- 
tional glitch size Afg/i/. Our results reinforce the bimodal 
distribution of the observed fractional gl itch sizes previ- 
ously reported by numerous authors (e.g., I Wang et al.ll2000l ; 
I Yuan et "all l2010al : lEspinoza et alj l2011al ). The first peak 
in this distribution lies around 2 x 10~ 9 and the second 
around 10 -6 . Our observations mainly contribute to the sec- 
ond peak. Because of our rather infrequent sampling, it is 
very difficult for us to detect glitches with Avg/v < 10 -9 . 



As noted bv lEspinoza et all (|2011al ) and others, the left edge 
of the distribution is strongly limited by observational selec- 
tion. The actual number of small glitches could be large and 
the lower peak in Figure[l2]may not even exist in the intrin- 
sic distribution. However the dip at Au s /v ~ 10~ 7 is clearly 
real and suggests that there may be two mechanisms that 
can induce a glitch event. As previously mentioned in it 
has been proposed that starquakes caused by the cracking 
of stellar crust may generate small glitches, whereas large 
glitches may result from the sudden transfer of angular mo- 
mentum from a crustal superfiuid to the rest of the star. 

The fractional glitch size is affected both by the size of 
the glitch and the pulse frequency of the pulsar. In the lower 
panel of Figure 1121 we plot a histogram of the frequency 
jump Af g . As shown in the figure, the distribution of Ai/ g 
also has a bimodal distribution or at least has a dip between 
the large and small glitches. It is interesting that the peak for 
large glitches is much narrower in Av g than it is in Ai/ g /V, 
whereas the converse is true for the lower-frequency peak. A 
large fraction of the peak at t he high end comes from just 
two pulsars, PSR J0537-6910 (|Middleditch et al.|[200^ ) and 
the Vela pulsar. These two pulsars have frequent glitches and 
most of them have Av g ~ 20 ^tHz. The pulse frequencies 
however differ by a factor of about six. 

Figure [13] shows the time sequence of glitch fractional 
sizes for seven pulsars where ten or more glitches were 
detected. This figure shows that the bimodal distribution 
of glitch sizes may be seen in individual pulsars as well. 
For example, most glitches in PSR J0537— 6910 and the 
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Figure 12. Histograms of the distributions of glitch fractional 
size Avg/v (upper panel) and of glitch size Au g (lower panel). 
The blank bars show the results from the previous work, while 
the shaded bars show the new detections from this work added on 
top of the previous results. For previously reported glitches, we 
use the ATNF glitch database which includes magnetar glitches. 



Vela pulsar (PSR B0833— 45) are large and similar in am- 
plitude, but much smaller glitches were occasionally seen. 
Although not quite so clear, similar behaviour is seen in 
PSRs J1341-6220, J1740-3015, J0631+1036 (many small 
glitches and only one large glitch) and J1801— 2304. It is ev- 
ident from Figure [T3]that we detected fewer small glitches in 
PSR J1341— 6220 compared to earlier results. This may be 
because the frequent occurrence of larger glitches obscured 
some small glitches in our data set which is not as well sam- 
pled as the earlier ones. 

In Figure 1141 we show a set of period-period-derivative 
(P-P) diagrams to present six quantities relevant to 
glitches. For previously published glitches, we refer to the 
ATNF Pulsar Catalogue glitch table. The six plotted quan- 
tities are a) the number of glitches detected in a given pulsar 
N g , b) the average number of glitches per year Ng, c) the 
fractional glitch size Avg/v, d) the glitch size Ai/g, e) the 
rms glitch fractional size and f) the rms fractional size nor- 
malised by the mean value for that pulsar. For sub-plots 
c) and d), if a pulsar has glitched more than once, then the 
largest value is plotted. In each P-P diagram, the size of the 
symbols (circle or triangle) is a linear function of the magni- 
tude of the given quantity; we adjusted the slopes and offsets 
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Figure 13. The time series of glitch events for the seven pulsars 
where ten or more glitches were detected. Data from the ATNF 
Pulsar Catalogue glitch table are indicated by circles, while those 
from this work are indicated by triangles. For each pulsar, its 
characteristic age is shown in parentheses. 

for the different functions to give appropriate sizes for the 
symbols. 

Table lA"3l gives the number of detected glitches Ng, the 
observation data span and the rate of glitches for known 
glitching pulsars Ng. The uncertainty of the glitch rate was 
estimated as the square-root of N s divided by the data span. 

We discuss these results in the following subsections. 

5.1.1 The number and rate of glitches 

Large numbers of glitches were observed in the pulsars with 
characteristic ages between 10 3 and 10 5 yr; the seven pul- 
sars that have been observed to show ten or more glitches 
are within this age region. If magnetic-dipole radiation is 
assumed, then the inferred surface magnetic field for the 
seven pulsars spans from ~10 12 to ~10 13 G. But there 
are some young pulsars that have not been observed to 
glitch, although they also have relatively long data spans. 
For instance, PSR J1513— 5908 has a characteristic age of 
^1.5 kyr; it is the youngest pulsar in our sample of 165 pul- 
sars. This pulsar has been observed for more than 28 years 
^Livingstone &: Kaspi|[201ll) with no evidence for any glitch 
activity. On the other hand, PSR J1119— 6127 has a similar 
characteristic age of ~1.6kyr and three glitches have been 
observed in its ~13-yr timing data span. 

The 23 glitches observed for PSR J0537-6910 oc- 
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Figure 14. P-P diagrams for glitch-related quantities of the a) number of detected glitches; b) average number of glitches per year; 
c) maximum fractional glitch size; d) maximum glitch size; e) rms fractional glitch size; and f) rms fractional size normalised by the 
mean. A circle indicates the parameter was obtained from the ATNF Pulsar Catalogue glitch table, whereas a triangle symbol indicates 
a parameter from this work. In the various plots, the seven pulsars exhibiting ten or more glitches are marked: 1 — PSR B0531+21 (Crab 
pulsar); 2 - PSR J0537-6910; 3 - PSR B0833-45 (Vela pulsar); 4 - PSR J1341-6220; 5 - PSR J1740-3015; 6 - PSR J0631+1036; 7 - 
PSR J1801-2304; and two magnetars: A - PSR J1048-5937 (IE 1048.1-5937) and B - PSR J1841-0456 (IE 1841-045). 
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curred in ~7.6yr, resulting in a large glitch rate of 
~3.0yr _1 ; this pulsar is the most frequently glitching pul- 
sar known. The second most frequently glitching pulsar is 
PSR J1740-3015. This pulsar has 32 glitches in ~25yr, 
giving a glitch rate of ~1.3yr - . As shown in sub-plot 
b) of Figure 1141 the positions of these two pulsars on 
the P-P diagram are quite different: the pulse period for 
PSR J0537-6910 is about a factor of 40 less than that of 
PSR J1740— 3015 and P is about an order of magnitude 
smaller. Hence the two pulsars are both relatively young, 
with characteristic ages of 5 kyr and 20 kyr respectively. 
The dipole magnetic field strength of PSR J1740-3015 
is more than an order of magnitude stronger than that 
of PSR J0537-6910. We found that PSRs J0729-1448, 
J1341-6220 and J0922+0638 also exhibit glitches more than 
once per year. For PSR J1341— 6220, the Parkes observations 
have detected 25 glitches during ~20yr. PSR J0729-1448 
has glitched four times during its data span of ~3.3yr. For 
PSR J0922+0638, one glitch has been d etected in a data 
span of just one year dShabanoval |2010T ). so the inferred 
high glitch rate is very uncertain. On the P-P diagram, the 
characteristc-age lines of 10 3 and 10 s yr together with the 
magnetic-field lines of 10 12 and 10 13 G define a region where 
pulsars are observed to exhibit more glitches and large glitch 
rates. 



5.1.2 Glitch sizes 

As shown in sub-plot c) in Figure [T4l large fractional glitch 
sizes Ai/g/v are generally observed in young pulsars with 
long periods (i.e., small u). The largest known, ~ 3.3 x 10~ 5 , 
was detec ted in PSR J1718— 3718 wh ich has a pulse period 
of 3.38 s dManchester fc Hobbsl I201H ) but a characteristic 
age of only 34 kyr and an inferred dipole field strength of 
7.4 x 10 13 G, one of the highest known for radio pulsars. 
Magnetars detected at X-ray wavelengths are also known to 
suffer large glitches and have long pulse periods and even 
higher inferred dipole fields. For example, PSR J1048-5937 
(IE 1048.1-5937) and PSR J1841-0456 (IE 1841-045) are 
both X-ray dete cted magnetars which have had glitches with 
Avjv > 10" 5 ijDib et al J 1 20081 l2009r i. 

Not surprisingly, the absolute frequency jumps Av & 
tend to be larger for the shorter-period (larger v) pulsars 
since they are not normalised by v. Sub- plot d) of Figure [14] 
shows the distribution of Au s values on the P-P diagram 
confirming this expectation. The largest frequency- jump was 
observed in the 24-ms pulsar J2022+3842; the pulse fre- 
quency of_t4Tis_s^urce_gamed_~78 /iHz in its MJD ~54675 
glitch jArzoumanian et all 120111 ). Sub-plot d) of Figure [Til 
also shows that the glitch size Au s is also correlated with 
characteristic age. Pulsars with young ages tend to have 
larger jumps, but this correlation breaks down for the very 
young pulsars such as the Crab pulsar (PSR B0531+21) and 
PSR J1513-5908. 



5.1.3 Glitch variability 

Sub-plot e) of Figure [14] shows the rms fluctuation in 
fractional glitch sizes. The magnetars PSR J1048— 5937 
(IE 1048.1-5937) and PSR J1841-0456 (IE 1841-045) 
show a wide range of glitch sizes, with Av g /v ranging from 



~ 1.4 x 10~ 6 up to ~ 1.6 x 10~ 5 . Although most of the 
glitches in the Vela pulsar are large, there have been two 
small glitches (see Figure I13[) and so the rms fluctuation of 
glitch size is relatively large. After the magnetars, large rms 
variations are found in the high-magnetic-field radio pulsar 
J1119-6127 and PSR J1838-0453. There seems a tendency 
for more variability in glitch size in older pulsars. This is also 
seen in sub-plot f) where the rms fluctuation has been nor- 
malised by the mean value to give an effective "modulation 
index" for glitch size fluctuations. However, this conclusion 
is not very certain as only a small number of glitches were 
detected in PSR J1838— 0453 and some other pulsars. Also, 
the results may be biased by the difficulty in detecting small 
glitches in noisy pulsars. In sub-plot f) there is more scatter 
in this plot though with the Crab pulsar being more promi- 
nent because of its relatively small mean glitch size. 



5.1.4 Implications for neutron-star physics 

In the angular-momentum-transfer model, a glitch is un- 
derstood as a sudden transfer of angular momentum from 
a more rapidly rotating interior superfluid to the neutron 
star crust fe.g.. lAlpar et al]|l98lf ). The angular momentum 
of a rotating superfluid is carried by an array of vortices; 
each vortex contains a quantised unit of angular momen- 
tum. In principle, if the rotation of a superfluid is slowing 
down, then the surface density of the vortices will be de- 
creasing, or, in other words, the vortices will migrate away 
from the rotation axis. However, in the solid crust of a neu- 
tron star, the vortices o f the neutron superfluid tend to pi n 
onto the nuclear lattice (|Alparl 19771 ; lEpstein fc Bavmf l988). 
The pinning between the vortices and the crystal lattice will 
not be broken until the force induced by the differential ro- 
tation (or the "Magnus" force) between the superfluid and 
the crust reaches a critical value, beyond which an avalanche 
process of unpinning and a substantial transfer of angular 
momentum could be triggered. 

Sub-plots c) and d) of Figure Q3] shows that large 
glitches, either relative or absolute, are mostly confined to 
the pulsars with characteristic ages smaller than 10 5 yr. The 
largest fractional glitch sizes are observed in long-period, 
low-i^ pulsars. This suggests that a larger fraction of the ex- 
cess angular momentum of the superfluid is transferred to 
the crust in these pulsars. For very young pulsars such as 
the Crab pulsar, glitch sizes are much smaller and somewhat 
less frequent. For a very young neutron star, as the result of 
the internal high temperature, the superfluid vortices could 
creep against the pinning energy barrier, preventing the for- 
mation o f pinning zones an d the sudden relea s e of p inned 
vortices (|Alpar et al.lll984al ). iRuderman et ail (|l998h have 
suggested that, because of tangling of the vortices and mag- 
netic flux tubes, as the neutron star spins down the core 
flux-tubes are pulled by the vortex lines, moving toward the 
stellar equator. The moving core flux-tube s lead to a build- 
up of the shear stress in the stellar crust JSrinivasan et all 
1990). Once the stress grows to exceed the yield strength of 
the crust, then the consequent cracking of the stellar surface 
may cause small peri od jumps in the rotation of the neutron 
star <|Rudermanll2009l ). This could be the mechanism for the 
observed small glitches in the Crab and other pulsars such 
as Vela and PSR J1740— 3015. It is also possible that such 
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Figure 15. Histogram for the recovery fraction Q of exponential 
decays. Data for the blank bars are from the literature as given 
in Table lA4l while the shaded bars show values from this work. 



processes may trigger the large-scale release of angular mo- 
mentum needed for glitches with large fractional glitch sizes. 



5.2 Glitch recoveries 

Pulsars show a variety of behaviours following a glitch. Part 
of the step change in both v and i> often recovers exponen- 
tially, in some cases with more than one identifiable time 
constant. Following this exponential recovery, a linear in- 
crease in v (decrease in slow-down rate \v\) is often observed. 
Normally this continues until the next glitch. Finally, appar- 
ently permanent changes m.v,v and/or v are sometimes left 
after the transient recoveries. We searched for evidence of 
these different types of recovery in all 107 observed glitches. 



5.2.1 Exponential recoveries 

In our sample, 27 glitches in 18 pulsars had an identifiable 
exponential recovery. The observed fractional recovery Q 
and recovery time constant Td for these glitches are given 
in Table IA2I All but two of the exponential recoveries are 
well modelled by a single exponential term. Two exponen- 
tial terms were required for glitch 3 in PSR J1119— 6127 
and glitch 1 in PSR J1803-2137. The largest Q ~ 0.84 
was detected in glitch 2 in PSR J1119-6127. Table lAHsum- 
marises the parameters for the previously reported exponen- 
tial recoveries. Multiple decays were observed for three pul- 
sars, the Crab pulsar, the Vela pulsar and PSR J2337+6151. 
However, the strong observational selection against observ- 
ing short-term recoveries in most pulsars needs to be recog- 
nisedH 

Figure [T5] shows the histogram of the fractional 
exponential recoveries Q. The observed values span a 
very wide range fr om the smallest ~ 0. 00014 observed in 
PSR J1841-0425 (|Yuan et all l2010al ) to the largest ~8.7 



6 Note that we do not consider t he 1.2-min recovery observed 
after the Vela MJD 51559 glitch bv lDodson et al.l <2002t) for two 
reasons: a) the parameters of the recovery are very uncertain and 
b) it is unique in that such short timescale recoveries cannot be 
observed for any other known pulsar. 



recently detected in th e young X-ray pulsar J1846— 0258 
^Livingstone et alj|2oioh . The histogram is clearly bimodal 
with a broad peak around 0.01 and another very close to 1.0. 
This strongly suggests that there are two different mecha- 
nisms for the exponential recovery following a glitch. 

To further explore the properties of the exponential re- 
coveries, all known values of Q and Td are shown in Figure 
1161 as a function of the fractional glitch size Avg/v, the pul- 
sar characteristic age and the surface magnetic-dipole field. 
In sub-plot a), it is striking that glitches with Q ~ 1 are 
found over the whole range of fractional glitch sizes, whereas 
those with small Q are only found in the larger glitches 
with Ais g /v > 10 -6 . Because of random period irregulari- 
ties, it will be difficult or impossible to detect exponential 
recoveries in many small glitches with Av g /u < 10~ 8 . How- 
ever, such decays should be detectable in most glitches with 
Aug/u ~ 10~ r . Therefore there appears to be a real absence 
of exponential recoveries with Q ~ 0.01 for glitches with 
Au g /v < 3 x 10~ 7 . It is also worth noting that, for many 
glitches, both large and small, no exponential recovery is de- 
tected (Table |A"2)) . so effectively Q ~ in these cases. These 
results suggest that the real physical distinction between 
pulsars with large glitches and those with small glitches (cf. 
Figure I12p is more complex with the glitch recovery param- 
eter also being important. 

Unlike Q, the distribution of the time constant Td rel- 
ative to Avgjv shown in sub-plot b) is more uniform with 
decay timescales from a few days to a few hundred days seen 
in glitches of all sizes. 

Sub-plot c) of Figure fTol shows that the large-Q glitch re- 
coveries are observed in both young and old pulsars, whereas 
the large glitches, which have low Q, are not seen in the 
very young pulsars. As discussed above, the large glitches 
primarily occur in pulsars with characteristic ages between 
10 4 and 10 s years. Sub-p l ot d) g ives some support to the 
suggestion of lYuan et al.l (|2010aT ) that there is a positive 
correlation of Td with pulsar characteristic age. There is 
a clear absence of short-term decays for the older pulsars 
that cannot be accounted for by observational selection. 
Sub-plots e) and f) show the dependence of Q and Td on 
surface dipole field strength. Since this is derived from P 
and P, as is the characteristic age, it is not surprising that 
these plots show basically the same dependences as sub- 
plots c) and d). However, they do highlight the fact that 
the m agnetars, with B B > 10 1 3 G, have high- Q glitch recov- 
eries l|Kaspi fc Gavriill 120031 ; iDib et al l l200Sl ; iGavriil et all 
l201ll ). The largest observed Q ~ 8.7, indicating a massive 
over-recovery, was detected in the young X-ray-detected but 
spin-powered pulsar PSR J1846— 0258, which has an inter- 
mediate period, ~0.326 s, a very high implied dipole mag- 
netic field, B B ~ 5 x 10 13 G and is located near the centr e of 
the supernova remnant Kes 75 (|Livingstone et al"1l2010l ). 



5.2.2 Linear recovery 

It has been well recognised that the long-term recovery from 
a glitch is generally dominated by a linear increase in z>, that 
is, a linear decrease in the slow-down rate or, equiva- 



lently, in P jDowndll98ll;lLvne et aL]ll996l : iLvne et alj|2000l : 
IWan'g et allhoOfj ; lYuan et al.ll2010ah . This effect is clearly 



seen in most of the v plots given in Figures 131- 111 1 especially 
for the larger glitches. The linear trend normally becomes 
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Figure 16. Dependences of the post-glitch exponential recovery parameters Q and on the observed glitch fractional size, the pulsar 
characteristic age r c and the dipole magnetic field B B : a) Q versus Ai>g/v, b) versus Av g /v, c) Q versus t c , d) Td versus r c , e) Q 
versus B B and r<i versus _B a . 



evident at the end of the exponential recovery and persists 
until the next glitch event. In cases where the exponential 
recovery is absent or has very low Q, e.g., PSR J1301— 6305 
(Figure [5]) the linear recovery begins immediately after the 
glitch. The rate of increase in v is quantified by fitting for 
v in the pre-, inter- and post-glitch intervals following the 
decay of any exponential recovery. Values of v obtained in 
this way are given in Table XK\\ and are plotted in Figure [T71 
for 32 pulsars. No significant it value was obtained for four 
glitching pulsars, in most cases because available data spans 
were too short. Of the 108 v values plotted, 11 are negative. 



These are for pulsars where timing noise is relatively strong 
and/or the available data spans are short. 

There is no doubt that these linear increases in i> are 
related to internal neutron-star dynamics and recovery from 
glitches. Firstly, the value of v often changes at the time 
of a glitch. Clear examples of slope changes at glitches are 
seen for the Vela pulsar (Figure [3J), PSR J1420-6048 (Fig- 
ure© and PSR J1709-4429 (Figure©. In other cases how- 
ever, e.g., PSR J1301-6305 (Figure© and PSR J1803-2137 
( Figure [IT))) , there is little or no slope change at a glitch. Sec- 
ondly, although a positive value of v is expected from normal 
magnetospheric braking, the observed values are generally 
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Figure 17. Linear decay rate v as a function of for 32 glitching 
pulsars. The symbols of a circle and triangle indicate positive 
and negative values, respectively. The solid line gives the v ey ± 
resulting from external or magnetospheric braking with n = 3 
and for v = 10 Hz. 



much larger. Pulsar braking is normally described by the 
braking index n, denned by 

where n = 3 for magnetic-dipole braking. This relation is 
shown in Figure [17] assuming n = 3 and for a typical young 
pulsar with v = 10 Hz. The magnetospheric or external con- 
tribution i> cxt is well below the observed values. Observed 
braking indices attributabl e to magnetospheric b raking are 
generally less than 3.0 (e.g.. lEspinoza et alj|2011bl ) which in- 
creases the discrepancy. However, for young long-period pul- 
sars, the discrepancy is less. For example, PSR J1119— 6127 
has a value of v for magnetic-dipole braking comparable to 
the minimum obser ved values (Table All) a nd its braking in- 
dex is close to 3.0 (jWeltevrede et aT. 20 111 ). However, there 
are significant changes in v at the glitches in this pulsar, 
with the largest observed value being more than a factor of 
two higher. So even in this case, it is clear that glitch-related 
phenomena contribute to the v value. 

Figure [TS] shows the change in slope of the linear recov- 
eries at glitches, Av, as a function of slow-down rate \v\ and 
glitch fractional size. Sometimes Ai> is solved for as part of 
the glitch fit (Table |A"2)) . but often the data spans used for 
this are too short to define Av well. Therefore, in order to 
enlarge the sample as much as possible, we took the differ- 
ence in v between each pair of the pre- and post-glitch solu- 
tions given in Table |A"T1 We found that, out of the available 
66 values, 35 are positive (53%) and 31 are negative (47%) 
in accordance with our expectation that the positive and 
negative changes would be approximately evenly balanced. 
Panel a) further illustrates that the changes are comparable 
in magnitude to the v values plotted in Figure [T71 implying 
that the amount of additional braking can change dramati- 
cally at a glitch. Panel b) shows that the amount of change 
is not strongly dependent on the glitch size, with both large 
and small glitches inducing a wide range of slope changes 
relative to the v value. Since large glitches usually occur 
in relatively young pulsars, most of which have a large \v\ 
(cf. Figure 1141 panel d) , it is not surprising that the largest 
values of AD occur in glitches with large Av s /v. 



5.2.3 Implications for neutron-star physics 

lAlpar et all (Il984al lbl. Il993f ) have suggested that the ob- 
served different types of glitch recoveries are a manifesta- 
tion of "vortex-capacitive" and "vortex-resistive" regions in 
the interior superfluid. For the former, the vortices decouple 
from the normal spin of the crust in a glitch but only share 
angular momentum at the time of the glitch, whereas for the 
latter, the vortices drift via a continuous pinning and unpin- 
ning, continuously coupling their angular momentum to the 
crust. At a glitch, the vortices in the resistive zones unpin 
and then gradually repin until the next glitch, resulting in 
the observed exponential and linear recoveries in v. The re- 
sistive region may contain several sub-layers, some of which 
may have a linear response and others a non-linear response 
to the glitch. In this context, a linear response means that 
the amplitude of the associated exponential recovery toward 
the steady state is proportional to the perturbation in the 
angular velocity of a particular sub-layer at the time of a 
glitch. This occurs when the internal termperature is high 
compared to the pinning energy and the differential angular 
velocity between t he superfluid layer and the crust is small 
jAlpar et alJll993h . On the other hand, in outer layers where 
the temperature is low compared to the pinning energy, the 
equilibrium lag is large and the response is not linear with 
the angular-velocity perturbation. Equilibrium is generally 
not reached before the occurence of the next glitch and there 
is an effectively linear incre as e in u . 

Following lAlpar et al] l| 19931 ). we label the linear- 
response layers that are responsible for the exponential re- 
coveries, i = 1,2,..., with moments of inertia 7i, and the 
nonlinear response layer, responsible for the linear recovery 
in layer A with moment of inertia 7a. For the exponential 
recoveries, the decaying part of Ai> is 

Ai/g 7i Ai/g 

Td,i I Td,i 

where 7 is the total moment of interia of the star. For the 
large glitches, the median observed value of Q is about 0.01 
(Figure fTrJK ) showing that the linear-response superfluid lay- 
ers contain about 1% of the total moment of inertia. 

The situation is clearly different for the pulsars for 
which Q ~ 1. In the simplest two-component models with 
a solid crust and a superfluid interior, where the super- 
fluid is weakly coupled to the crust, Q « 7 n /7, where 7 n 
is th e total moment o f inertia of the superfluid neutrons 
(e.g.. lBavm et af]|l969l ). So one explanation for the high-Q 
events is that a large fraction of the total stellar moment of 
inertia is in the form of superfluid neutrons that are weakly 
coupled to the crust. The degree of coupling is evidently 
quite variable, with decay time constants ranging from a 
few days for the Crab pulsar to tens or hundreds of days 
for other pulsars. This explanation for high-Q decays can 
apply regardless of the mechanism for the glitch itself; this 
may be related to the fact that these decays are observed to 
follow both large and small glitches. The observed overshoot 
in the v ery young pulsar J1846— 0258, with an apparent Q 
of ~8.7 (|Livingstone et aill2010l) . is clearly anomalous with 
no other similar oversho ots being observed. As suggested by 
iLivingstone et all {2010.:) it is likely that this behaviour re- 
sulted from a glitch-induced change in the external torque. 

In the Alpar et al. models, the observed approximately 
linear recoveries in v are related to the properties of an outer 
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Figure 18. The change in slope of linear i> recoveries at glitches, Ai>, as a function of a) \0\ and b) glitch fractional size Au g /u. Circles 
and triangles indicate the postive and negative values, respectively. 
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Figure 19. The observed glitch-related i>i nt as a function of the 
ratio of | i> | and the mean interglitch interval (r g ). The solid line 
shows a fitted straight line which has a slope very close to 1.0 (see 
text). 



superfiuid layer in which the coupling is very weak so that 
the rotational lag betwe en the crust and the s uperfiuid is 
very large. For this case. lAlpar fc Bavkal (2006) show that 
the Ant related to glitch recovery is given by 



fint 



(6) 



where r g is the interglitch interval and I a is the inertial 
moment of layer A. Figure [191 shows Ant = A>bs — iv 2 /v ver- 
sus |^|/{7" g ), where (T g ) is the mean interglitch interval (Ta- 
bic [X3j. A correlation is clearly seen, with an unweighted- 
least-squares fit giving 



A nt = 10- 2 - 8 ( 1 - 4 '(|/>|/(r g )) 1 - °( 7 ). 



(7) 



Remarkably, the slope of the fitted line is equal to the ex- 
pected value of 1.0 based on Equation Furthermore, the 
proportionality constant, ~ 10~ 2 ' 8 ~ 0.0016 is comparable 
to the median value of Au g /i>, 0.0034 (Table |A"2)). although 
this quantity has a large scatter. This is consistent with the 
non-linear respons e models where both of these quantitie s 
are equal to I A /I jAlpar et al.lll993l ; lAlpar fc Bavkaf 2006). 




o 



Days from MJD 52500.0 



Figure 20. Timing residuals for PSR J1539-5626 after fitting 
up to v are shown in the top panel. Points marked with * are 
from digital filter bank observations. The bottom panel shows the 
variations i>. 



5.3 Slow glitches 

So-called "slow glitches" have been observed in a number 
of pulsars. In these events, v rises sharply and declines over 
the next few hundred days roughly to its pre-glitch state. 
This results in a persistent increase in v relative to the pre- 
glitch variation. Slow glitches were first observed in PSR 
J182 5 -0935 (PSR B182 2-09) l|Zou et allbool ; IShabanoval 
120071 : lYuan" et al.ll2010ah and have bee n reported in several 
other pulsars by lYuan et all l|2010al ). iHobbs et aD (|2010h 
suggested that these slow glitches are a manifestation of the 
discrete states in spi n -down rate first se en in PSR B1931+24 
ijKramer et al.ll2006h . rLvne et all ()2010h further showed that 
the discrete spin-down states were correlated with pulse 
shape changes, implying that the slow-glitch glitch phe- 
nomenon has a different origin to normal glitches for which 
the pulse shape changes are not expected. Never-the-less 
they are discrete events which result in a step change in 
spin frequency, so it is reasonable to label them as "slow 
glitches" . 
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Our data for PSR J1825-0935 (Figure [TTJ do not 
clearly show the slow glitches since there were insuf- 
ficient observations at the relevant times. However, for 
PSR J1539— 5626, we see slow-glitch features as shown in 
Figure [20l The top panel shows the timing residuals for this 
pulsar showin g the quasi-sinusoid al features as observed in 
many pulsars l|Hobbs et al.|[201ol ). The bottom panel shows 
the episodic increases in v that characterise the slow-glitch 
phenomenon. The two states have v ~ —8.19 x 10 -14 s -2 
and v ~ —8.15 x 1 0~ 14 s~ 2 , a relative variation of just 0.5% 
m v. iLvne et al l (|2010l ) found values of Avjv of between 
0.3% and 45% in different pulsars. Furthermore, they found 
a relationship between the size of the change in v and the 
size of the pulse-shape change between the two states. The 
small percentage change in v for PSR J1539— 5626 therefore 
suggests that any correlated pulse shape changes will also 
be small and so far we have been unable to find convincing 
evidence for the changes in either pulse shape or polarisation 
propertities. 



6 CONCLUSION 

In this paper, we reported the results of a search for glitch 
events in the timing residuals of 165 pulsars, covering a to- 
tal data span of 1911 yr. A total of 107 glitches was de- 
tected in 36 pulsars of which 13 were not previously known 
to glitch and 46 glitches are new discoveries. As constrained 
by our observational sampling, glitches with Au g /u < 10~ 9 
are difficult for us to detect, and 22 events that have pre- 
viously been reported were missed. Similar difficulties also 
occur for the detection of post-glitch exponential recoveries 
with timescales < 20 d. However, our observations do re- 
veal exponential recoveries mostly with timescales of a few 
tens of days for 27 glitches in 18 pulsars. A linear increase 
in z> is clearly observed following most glitches. Linear in- 
creases, presumably related to a previous unseen glitch, were 
also seen before the first observed glitch in most pulsars. To 
quantify the linear recoveries as accurately as possible, the 
solutions for v (see Table [AT]) were measured after any ob- 
served exponential recoveries were essentially complete. Of 
the 108 v measurements obtained, 97 are positive; the 11 
negative ones are generally from short data spans and/or 
for pulsars with strong timing noise. 

With the contribution of 46 new glitches, the observed 
bimodal distribution of the glitch fractional size has been re- 
inforced, implying that there are two different glitch mech- 
anisms, possibly the starquake and the vortex pinning - 
unpinning theories. Post-glitch exponential recoveries have 
been observed over a wide range of fractional glitch size and 
pulsar age. Large recovery fractions Q have been seen in 
small glitches in both young (e.g. the Crab) and old (e.g. 
PSR B0525+21) pulsars and also in large glitches in young 
pulsars (e.g. PSR J1119— 6127) and in magnetars. Small val- 
ues of Q, typically around 0.01, are more commonly ob- 
served in large glitches. A bimodal distribution shown in 
the histogram for Q has also clearly been seen. Moreover, 
decay timescales have been observed to show some pos- 
itive correlation with pulsar age. Figure [TT] shows that the 
inter-glitch v has a strong correlation with the slow-down 
rate \v\ and is generally much larger than the expectation 
from a magnetic-dipole braking. The excess decay in brak- 



ing is clearly related to glitches and is consistent with the 
predictions of a model based on the properties of a weakly 
coupled superfluid in the outer layers of the neutron star 
I Alpar fc Bavkalll2006l '). 

It is very clear that the true distribution of the glitch 
fractional size at the low end is not well determined (cf. 
Figure [T2"|) . Glitches with Av g /is < 10 -9 are at the lower 
limit of the detectability for most timing programs. The ob- 
servational sampling, timing accuracy and intrinsic timing 
noise also hamper the detection of exponential recoveries 
with very short timescales, leaving an incomplete sample for 
short time constants. Further studies on these issues require 
intensive timing observations for the glitching pulsar pop- 
ulation supported by simulations to better reveal the true 
distributions of glitch-related phenomena. 
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APPENDIX A: TIMING SOLUTIONS AND 
GLITCH PARAMETERS 

Timing solutions incorporating pre- and post-glitch pulse 
parameters and glitch parameters, obtained from Tempo2 
fits for the 36 glitching pulsars analysed in this paper are 
provided here in Tables |A"T1 and [A"2l 

To support the discussion in ij5] Table IA3I summarises 
the number of observed glitches, iV g , the observing range 
and the derived glitching rates N s for the known glitch- 
ing pulsars. Table 1X41 gives previously reported exponential- 
recovery parameters. 

ASCII machine-readable versions of the four tables in 
this Appendix are provided as an on-line supplement. 
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Table Al. Pre- and post-glitch timing solutions for 36 glitching pulsars. 
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- 


49425 


49364 




49488 


6 





.25 


5.64[3] 






1- 


2 


88776198459(6) 


-0 


.20845(3) 


- 


49649 


49560 




49739 


9 





.26 


2.36[6] 






-2 


2 


88771636557(6) 


-0 


.2071686(11) 


-2.02(12) 


52195 


51526 


— 


52865 


21 





.95 


355[17] 






2- 


2 


88768278363(6) 


-0 


.2073986(12) 


0.38(7) 


54071 


52998 




55145 


65 


3 


.70 


2790[61] 


.11016 


-5857 


-1 


9 


3126315989(4) 


-6 


.990896(13) 


78(1) 


51913 


51299 




52527 


97 


3 


.28 


1830[93] 






1-2 


9 


3115113267(6) 


-6 


.994765(11) 


122.4(5) 


53791 


52571 


— 


55011 


134 


15 


.3 


27700[130] 






2- 


9 


3106478766(6) 


-7 


.00431(5) 


133(18) 


55250 


55072 




55429 


19 





.78 


64.8[15] 


J1048 


-5832 


-1 


8 


0873027360(12) 


-6 


.27372(4) 


76(5) 


48418 


47910 




•1 S 9 2 8 


60 


10 


.6 


52500[56] 






1-2 


8 


08699243910(13) 


-6 


.27060(14) 


- 


48991 


48957 




49025 


9 





.09 


3.85[6] 






2-3 


8 


0863750399(8) 


-6 


.27954(3) 


98(3) 


50172 


49559 


— 


50786 


61 


7 


.22 


148000[57] 






3-4 


8 


0854473241(13) 


-6 


.26780(3) 


122(3) 


51894 


51093 




52696 


61 


16 


.4 


1000000[57] 










fl -ZA T A A '3 OTP. < 1 1 1 






181(7) 




52771 




53653 




6 


.29 


96600 [3 1] 






5-6 


8 


0842772583(18) 


-6 


.26424(9) 


223(12) 


54082 


53680 




54486 


34 


7 


.37 


393000(30] 






6- 


8 


0838894263(12) 


-6 


.28529(6) 


98(12) 


54843 


54505 


— 


55183 


34 


4 


.42 


214000(30] 


J 1052 


— 5954 






5372526159(4) 




6134 1(11) 






54220 




54485 




1 


.25 


8 34[9] 






1- 


5 


53721546844(16) 


-0 


.617277(7) 


35(2) 


55096 


54733 




55461 


27 





.98 


3.37[23] 


,11105 


-6107 


-1 


15 


8248209436(12) 


-3 


.95900(6) 


-108(8) 


49995 


49589 




50402 


41 


2 


.54 


2590[37] 






1-2 


15 


8245012491(16) 


-3 


.96232(4) 


38(6) 


50942 


50434 




51451 


108 


6 


19 


17400[104] 






2-3 


15 


823737653(2) 


-3 


.96361(4) 


13(2) 


53217 


51744 




54690 


99 


28 


.0 


661000(95] 






3-4 


15 


8231268298(20) 


-3 


.96587(13) 


-11(30) 


55002 


54733 




55272 


38 


3 


.38 


13200(34] 






4- 


15 


8230116650(6) 


-3 


.96652(9) 


319(81) 


55382 


55304 




55461 


11 





.13 


35.9(7] 


J1112 


-6103 


-1 


15 


3936495872(5) 


-7 


.45479(4) 


-24(11) 


51055 


50850 




51261 


50 





.67 


6.61(46] 






1-2 


15 


3927981288(7) 


-7 


.471152(13) 


222(1) 


52417 


51529 




53307 


58 


5 


.26 


904(54] 






2- 


15 


.3913344780(9) 


-7 


.46853(4) 


144(5) 


54712 


54220 




55206 


42 


3 


.21 


344(38] 


J1119 


-6127 


-1 


2 


45326884390(8) 


-24 


.211160(7) 


700(2) 


51121 


50852 




51392 


57 


1 


.08 


7.46[53] 






1-2 


2 


4507184171(3) 


-24 


.142246(5) 


626.2(4) 


52342 


51405 




53279 


149 


21 


.04 


3800(145] 






2-3 


2 


447638860(2) 


-24 


.0540(1) 


814(14) 


53821 


53423 




54220 


33 


19 


.5 


3230(29] 






3- 


2 


444874218(3) 


-23 


.95777(11) 


773(17) 


55154 


54733 






56 


64 


.8 


11000(52] 


J1301 


-6305 


-1 


5 


4190840865(3) 


-7 


.82942(1) 


248(2) 


51420 


50941 




51901 


62 


3 


.61 


17.1[58] 










4182027029(12) 




S/i r\oo < a \ 






52145 




53371 




17 


.3 


4fc O'-i: [OOJ 






2- 


5 


417207268(1) 


-7 


833490(18) 


276(2) 


54249 


53395 


— 


55104 


64 


20 


.2 


580(60] 


J1341 


— 6220 






1 TOO ipOPOO 1 K \ 




TP. TP. ^ < P. \ 






49540 




49763 







.47 








1-2 


5 


1728273829(3) 


-6 


.76580(19) 




49837 


49787 




49888 


6 





.30 


11.3(3] 






2-3 


5 


1727579534(9) 


-6 


.7529(14) 


- 


49956 


49920 




49992 


7 





.65 


8. 79 [4] 






3-4 


5 


17263893435(17) 


-6 


.77137(3) 


167(9) 


50174 


50026 




50322 


16 





.36 


4.75(12] 






4-5 


5 


.1724951732(7) 


-6 


.76951(12) 


356(93) 


50420 


50341 




50501 


13 





.42 


9.74[9] 






5-6 


5 


1723882419(4) 


-6 


.76865(7) 


484(74) 


50603 


50537 




50670 


17 





.31 


8.02(13] 






6-7 


5 


1721912539(7) 


-6 


.76959(5) 


207(20) 


50946 


50760 




51133 


30 


1 


.29 


158(26] 






7-8 


5 


1719833589(3) 


-6 


.76796(4) 


308(13) 


51303 


51155 


— 


51451 


8 





.41 


10.9[4] 






8-9 


5 


1716336852(6) 


-6 


.76774(8) 


-63(36) 


51911 


51782 




52041 


16 





.93 


68.5(12] 










1714729568(9) 










52145 




52235 







.48 








10-11 


5 


171283363(4) 


-6 


.7866(3) 


-803(63) 


52518 


52266 




52771 


17 


11 


.8 


11700(13] 






11-12 


5 


170994398(12) 


-6 


.7534(8) 


-2120(269) 


53013 


52804 




53222 


15 


32 


.0 


188000(11] 






12-13 


5 


1708006160(7) 


-6 


.7484(3) 


- 


53347 


53241 




53455 


9 


2 


.01 


640(6] 






13-14 


5 


170541851(19) 


-6 


.7618(8) 


-177(181) 


53799 








19 






26.4(11] 






14-15 


5 


1702517491(3) 


-6 


.76414(4) 


365(14) 


54297 


54144 




54451 


15 





.58 


189(11] 






15-16 


5 


1700342469(7) 


-6 


.76214(7) 


287(22) 


54672 


54486 




54860 


15 


1 


.71 


189(11] 






16-17 


5 


169858206(3) 


-6 


.7592(4) 


3044(251) 


54976 


54881 




55072 


10 


1 


.90 


362(6] 






17- 


5 


16968689489(14) 


-6 


.764590(16) 


109(5) 


55282 


55104 




55461 


20 





.45 


13.5(16] 


J1412 


-6145 


-1 


3 


17232313958(4) 


-0 


.9928240(13) 


3.80(17) 


51353 


50850 




51858 


45 





.79 


2.59(41] 






1- 


3 


17212789172(5) 


-0 


.9969084(6) 


6.06(3) 


53887 


52314 




55461 


101 


4 


.73 


110(97] 


.11413 


-6141 


-1 


3 


5012521433(5) 


-4 


.08114(8) 




51055 


50850 




51261 


26 


7 


.20 


36.7(23] 






1-2 


3 


5011397601(19) 


-4 


.0842(8) 




51374 


51294 




51454 


9 


5 


.52 


21.8[6] 






2-3 


3 


50105387384(14) 


-4 


.08410(5) 




51627 


51472 




51782 


24 





.90 


3.02(21] 






3-4 


3 


5009433012(4) 


-4 


.08332(7) 


-350(40) 


51941 


51844 




52039 


16 





.79 


3.21(12] 






4-5 


3 


5007434868(4) 


-4 


.084389(14) 


107(3) 


52515 


52145 




52886 


27 


3 


.05 


60.9[23] 






5-6 


3 


5005683237(18) 


-4 


.0871(3) 


-1648(218) 


53012 


52925 




53101 


6 


1 


.61 


35.9[2] 






6-7 


3 


5003276299(11) 


-4 


.08149(4) 


53(4) 


53708 


53150 




54268 


36 


17 


.5 


1470(32] 






7- 


3 


4999222810(7) 


-4 


.077751(19) 


33(2) 


54882 


54303 




55461 


54 


14 


.0 


630(50] 


.11420 


-6048 


-1 


14 


6675125451(3) 


-17 


.83946(3) 


651(7) 


51311 


51100 




51523 


26 





.37 


15.5(22] 






1-2 


14 


666147261(2) 


-17 


.83559(11) 


964(12) 


52207 


51678 




52737 


80 


13 


9 


17000(76] 






2-3 


14 


6644360883(7) 


-17 


.85639(4) 


533(6) 


53336 


52957 




53716 


58 


2 


.51 


618(54] 






3-4 


14 


6631475381(9) 


-17 


.85577(3) 


945(5) 


54183 


53734 




54633 


38 


2 


.20 


599(34] 






4-5 


14 


6618531329(16) 


-17 


.84922(8) 


1552(14) 


55031 


54672 




55391 


46 


4 


.95 


2.91(42] 






5- 


14 


6612348736(20) 


-17 


.90(12) 




55445 


55429 




55461 


5 





.15 


4.40(2] 



Glitches in 36 southern pulsars 



Table Al. — continued 



PSR J 




(s" 1 ) 


(10~ 12 s~ 2 ) 


(io- 24 s" 3 ) 


Epoch. 
(M.ID) 


Data range 
(M.JD) 


ToAs 


Rms res. 
(ms) 


x r [d.o.f] 


T1452 


6036 


-1 


6.451956870952(9) 


-0.060343(1) 


- 


54635 


54220 


55051 


35 


0.16 


52.2[32] 






1- 


6.451953766111(19) 


-0.060447(5) 


- 


55266 


55072 — 


55461 


17 


0.12 


33.4(14] 






-1 


5.571461007545(8) 


-0.08520253(5) 


0.0659(19) 


52608 


OU D d y — 




117 


0.29 


643[113] 






1 — 


5.571444247574(7) 


— 0.0852032(9) 




54885 


54566 — 


55206 


28 


0.09 


237(25] 






-1 


11.87624700867(4) 


-1.938066(6) 




51448 






31 


0.15 


4.98[28] 






1 — 


11.8758384521(3) 


— 1.946724(4) 


14.85(15) 


54060 


52659 — 


55461 


99 


5.11 


6700[95] 


J1614- 


5048 


-1 


4.31760554(1) 


-9.1754(2) 


366(13) 


48848 




49 784 


100 


447 


2740000(96] 






1 — 2 


4.316214218(3) 


— 9.22347(7) 


118(5) 


50635 






171 


95.5 


632000(167] 








4.314824505(11) 


-9.1601(3) 


366(32) 


52385 


51782 


52989 


43 


159 


3800000(39] 






2- 


4.313373849(7) 


-9.17841(7) 


-48(5) 


54248 


53036 — 


55461 


100 


282 


4330000(96] 


J1646- 


4346 


-1 


4.3173283630(18) 


-2.095392(11) 


32.8(3) 


50857 


47913 — 


53803 


258 


294 


1080000(254] 






1- 


4.3166537681(7) 


-2.091929(20) 


109(2) 


54610 






47 


13.1 


718[43] 


J1702- 


4310 


-1 


4.15729055112(16) 


-3.864815(4) 


40.25(14) 


52498 


51223 — 


53774 


76 


7.49 


1100[72] 






1- 


4.15649409771(9) 


-3.868451(3) 


51.1(4) 


54941 


54421 


55461 


36 


1.15 


10.1[32] 






-1 


9.76127199222(13) 


-8.863764(4) 


123.9(8) 


48327 


47910 


48746 


46 


0.56 


487[42] 






1-2 


9.7597767546(8) 


-8.853391(18) 


179.7(7) 


50305 


49160 




168 


17.5 


385000[164] 






2-3 


9.758422208(3) 


-8.8516(1) 


595(11) 


52091 


51524 


52659 


36 


17.1 


1280000(32] 






3-4 


9.7570513620(17) 


-8.85369(4) 


244(4) 


53919 


53150 — 


54689 


60 


18.7 


1120000(56] 






A — 


9.7560834373(5) 


-8.86295(3) 


253(7) 


55219 


54932 — 


55507 


19 


1.21 


3890[15] 




3825 


— 1 


13.39192401122(9) 


— 2.3679369(7) 


25.38(3) 


52890 


50878 


54903 


146 


2.85 


1500[142] 






1- 


13.39144806554(9) 


-2.362796(5) 


16(2) 


55219 






18 


0.14 


7.09[14] 




3350 


-1 


7.170087711(1) 


-4.35574(2) 


85(2) 


51303 






69 


11.5 


35900(65] 






1 — 


7. 1690354505(5) 


-4.361361(6) 


50.5(3) 


54155 


52805 


^ r r 07 


90 


14.2 


36200(86] 


J1731 


4744 


-1 


1.20511216181(5) 


-0.2382630(12) 


18.89(12) 


48773 


48184 


49363 


21 


0.41 


8.58(17] 






1-2 


1.20508590674(6) 


-0.2376701(20) 


2.09(17) 


50059 


49415 — 


50703 


47 


3.18 


6330(43] 






2-3 


1.20505448010(3) 


-0.2375479(7) 


1.40(4) 


51590 


50722 — 


52458 


60 


1.53 


2230(56] 






3-4 


1.20502077952(3) 


-0.2375618(16) 


3.9(4) 


53239 






24 


0.50 


108(20] 






4 — 


1.204993924876(15) 


— 0.2374356(6) 




54548 






64 


1.96 


23500(61] 


_ 


3137 


-1 


2.21990926991(9) 


-0.68317(5) 


- 


54286 


54221 


54351 


10 


0.33 


0.75(7] 






1- 


2.21987416781(9) 


-0.684329(8) 


- 


54930 


54353 


55507 


49 


8.16 


141(46] 


J 1740- 


3015 


— 1 


1.648184042812(13) 


- 1 .265524(4) 


- 


50798 


50669 — 


50927 


42 


0.26 


82.3(39] 






1-2 


1.64809155129(19) 


-1.266547(13) 




51665 






27 


19.8 


399000(24] 






2 — 3 


1.6479812921(9) 


— 1.26564(5) 


125(98) 


52675 


52361 


52989 


19 


13.0 


179000(15] 






3-4 


1.64786917225(7) 


-1.266266(2) 


13.49(18) 


53728 


53036 — 


54420 


48 


2.16 


26600(44] 






4-5 


1.64774895802(16) 


-1.264817(7) 


11(2) 


54828 


54450 — 


55206 


37 


3.65 


134000(33] 






5- 


1.64769409872(4) 


-1.266147(5) 


70(3) 


55370 






17 


0.24 


625(13] 


J1801- 


2304 


-1 


2.4051703923(3) 


-0.65365(4) 


- 


48176 


47911 — 


48442 


16 


3.19 


23.0(13] 






1-2 


2.40512164856(9) 


-0.653541(3) 


2.8(3) 


49054 






61 


3.26 


21.5(57] 






2-3 


2.40507528156(11) 


-0.65348(5) 




49878 


49730 — 


50026 


11 


0.93 


5.07(8] 






3-4 


2.4050548970(3) 


-0.65345(4) 


29(19) 


50240 


50117 — 


50363 


25 


1.19 


7.83(21] 






4-5 


2.40502945211(8) 


-0.653426(7) 


40(2) 


50694 






83 


1.32 


8.79(79] 






5 — 6 


2.4049822461(6) 


— 0.65254(3) 


12(3) 


51531 


51021 


52041 


20 


7.93 


464(16] 






6-7 


2.40491876991(6) 


-0.653179(5) 




52684 


52145 — 


53223 


17 


2.22 


45.8(14] 






7-8 


2.40482774696(6) 


-0.653049(1) 


2.44(7) 


54318 


53307 — 


55330 


73 


3.72 


37.5(69] 






8- 


2.40476474326(9) 


-0.65292(5) 


" 


55435 






5 


0.22 


0.24(2] 


J1801- 


2451 


-1 


8.0071549008(15) 


-8.17840(11) 


369(32) 


49171 


48957 — 


49386 


17 


2.71 


1020(13] 






1-2 


8.006538717(1) 


-8.19043(3) 


387(3) 


50064 




5(1646 


81 


8.59 


8550(77] 






2 — 3 


8.005640187(4) 


— 8.18049(8) 


505(9) 


51348 


OUDOD 




97 


41.5 


234000(93] 






3-4 


8.004687982(3) 


— 8.19041(14) 


47(36) 


52737 


52484 — 


52990 


19 


4.77 


4330(15] 






4-5 


8.0039126298(19) 


-8.17244(5) 


146(4) 


53834 


53036 — 


54633 


58 


26.5 


73400(54] 






5- 


8.0029848506(6) 


-8.19035(3) 


266(7) 


55182 


54857 — 


55507 


35 


2.16 


343(31] 


J1803- 


2137 


-1 


7.48329875057(11) 


-7.48842(10) 




50709 


50669 — 


50750 


16 


0.12 


9.17(13] 






1-2 


7.4820940765(8) 


-7.4904029(17) 


244(2) 


52602 


51782 — 


53423 


58 


11.8 


101000(54] 






2- 


7.4807378694(16) 


-7.50335(4) 


288(4) 


54739 


53949 — 


55530 


58 


25.4 


192000(54] 


J1809- 


1917 


-1 


12.08488375958(18) 


-3.727797(2) 


37.74(12) 


52012 


50782 — 


53242 


57 


1.99 


365(53] 






1- 


12.08405907802(11) 


-3.7299474(3) 


33.63(16) 


54632 


53734 — 


55530 


62 


1.04 


62.6(58] 


J1825- 


0935 


-1 


1.3003898596(5) 


-0.088534(7) 


-1.4(5) 


52978 


51844 — 


54112 


69 


46.3 


578000(65] 






1- 


1.30037753992(5) 


-0.0888195(16) 


3.60(18) 


54608 


54144 — 


55073 


16 
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Glitch epoch determined by ph ase fit in this work. Refe renc es are given for glitch ep ochs adopted from previousl y pu blished work. Oth er values are determine d fro m our data sets. 
Refere nces: 1 - iFlanagar] l|l996l l: 2 - iDodson et alj lj2002h ; 3 - iDodson et al.ll|2004h : 4 - iFlanagan fc Buchnen (|2006h : 5 - lUramal (|2002l l: 6 - fEspinoza et"aD l|2011al ): 7 - iKrawczvk et alj 
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Table A3. Number of observed glitches and their mean rate for known glitching pulsars. 



Data span 
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2 
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Table A3. — continued 



Data span 
(MJD) 



References 
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1 


3(3) 
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Table A4. Previously reported parameters for the exponential recoveries that are not covered by our observations. 



PSR J 
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Age 




Gl. Epoch 






Q 
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(kyr) 


(10 ^ G) 


(MJD) 
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4U 01424-61 
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